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This  report  is  a  summary  of  a  theoretical  investigation  of  three  ways  to  avoid 
the  severe  problems  caused  by  recoil  in  tne  implanting  of  Propellant  Emplaced 
Anchors  (PEA's).  It  was  specifically  requested  that  a  launching  system 
modeled  after  a  recoilless  rifle  be  evaluated..  While  doing  this,  the 
conclusion  was  reached  that  the  recoilless  rifle  approach  was  more  complicated 
than  necessary,, and  that  the  desired  performance  could  be  achieved  more  easily 
by  a  direct  rocket.  The  sponsors  had  suggested  a  water  Jet  approach  using  cold 
high  pressure  gas  to  supply  the  energy,  instead  of  propellant  gas.  This  turned 
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out  to  be  impractical  using  the  compressed  air;  but  it  became  very  practical 
using  propellant  gases. 

All  of  the  systems  were  analyzed  and  simulated  by  computer  models.  All 
three  systems  appear  reasonable  and  they  all  solve  the  recoil  problem.  However, 
the  recoil! ess  and  the  direct  rocket  have  a  common  problem.  Their  require¬ 
ment  of  a  rocket  motor  with  very  high  thrust  (1.6  million  pounds)  and  a  very 
short  burning  time  (40  milliseconds)  has  no  parallel  in  current  rocket 
technology.  Apparent  solutions  to  this  problem  would  require  very  large 
structures.  In  any  case,  their  development  promises  to  be  expensive  and 
lengthy.  Fortunately,  the  water  jet  does  not  share  this  problem.  Its 
development  should  be  timely,  and  costs  should  be  reasonable. 
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FOREWORD 


The  work  on  propellant  emplaced  anchors  (PEA)  reported  here  was  funded  by 
the  Naval  Civil  Engineering  Laboratory  (NCEL)  under  Amphibious  and  Advanced 
Base  Technology,  PE  62760N,  Subproject  YF  60.536.  The  purpose  of  the  work 
covered  b_/  this  report  is  to  analyze  and  evaluate  three  candidate  propulsion 
systems  for  PEA's  to  replace  the  high  pressure  gun  propul sor  now  used.  These 
alternate  propulsion  systems  all  have  the  capability  to  propel  the  PEA's  and 
to  eliminate  the  recoil  problems  that  trouble  existing  PEA's.  All  of  the 
propulsion  systems  are  recoil  less.  All  three  propulsion  systems  were  analyzed 
in  depth,  and* the  analyses  were  implemented  into  conpiter  programs  that  can 
simulate  thei n, performance.  The  results  from  these  computer  programs,  in 
turn,  were  used,  to  make  estimates  of  the  size  and  weights  of  tne  competing 
systems.  The  fi\pal  selection  of  the' water  jet  propulsion  system  was  based  on 
size,  weight,  ane)  anticipated  short  development  time. 

\  Approved  by: 

\  . 

\  J.  F.  PROCTOR,  Head 

•  ■  •  \  Energetic  Materials  Division 
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A.  INTRODUCTION 

I.  THE  PROBLEM:  Newton's  Third  Law  has  teen  formulated  in  many  ways, 

but  they  all  have  the  same  meaning--all  forces  are  produced  in  equal  and 

opposite  pairs.  The  action  and  reaction  are  exactly  equal  in  magnitude  and 

exactly  opposite  in  direction.  Anyone  who  has  fired  a  gun  or  rifle  much 

larger  than  .22  caliber  understands  Newton's  Third  Law.  The  high  pressure 

gases  from  the  burning  propellant  act  both  on  the  base  of  the  bullet  and  the 

base  of  the  gun  breech,  producing  equal  forces  on  both.  Since  the  bullet  is 

free  to  move,  the  force  drives  it  down  the  barrel  and  out  the  muzzle.  The 

person  firing  the  gun  is  not  usually  free  to  move  very  far  because  he  is 

coupled  to  the  earth  (Figure  1).  The  recoil  force  is  transmitted  through  the 

shooter's  body  to  the  earth.  If  the  bullet  is  a  .30  caliber  with  a  weight  of 

150  grains  and  a  muzzle  velocity  of  3,000  ft/s,  the  recoil  force  will  change 

-24 

the  earth's  velocity  by  4.7  x  10  ft/s.  This  will  move  the  earth  one  inch 
in  five  trillion  centuries.  (The  reader  should  not  worry  about  the  accum¬ 
ulated  effects  of  shooting  on, the  earth's  trajectory,  because  air  friction  and 
the  final  impact  produce  equal  and  opposite  forces  which  exactly  counteract 
the  velocity  increment  of  the  earth.)  If  the  shooter  is  suspended  free  of  the 
earth,  his  much  smaller  mass  will  be  accel  jted  to  a  speed  of  four  inches  per 
second. 

The  previous  paragraph  discussed  the  'Cts  of  a  man  firing  a  .30  cal¬ 
iber  rifle  while  standing  oh  the  earth  and  while  suspended  in  the  air.  The 
effects  are  small  but  different.  Larger  caliber  guns  require  much  larger  and 
stronger  supports  than  a  rifleman  to  keep  the  gun  coupled  to  the  earth.  Our 
problem  is  greater  yet.  Propellant  Emplaced  Anchors  (PEAsj  comprise  a  system 
of  large  anchors  for  semipermanent  mooring  of  ships.  The  anchors  are  shot 
into  the  sea  floor  to  increase  their  effective  lifetimes'  and  their  load  capa- . 
cities.  There  are  four  load  ranges  from  10,000  pounds  to  300,009  pounds. 
Depending  on  their  holding  capacities  and  the. type  of  sea  floor  they  are 
intended  for,  the  anchors'  weights  range  from  160  pounds  to  6,890  pounds,  and 
their  required  velocities  range  from  360  to.  520  ft/s.  These  are  massive 
"bullets."  Projecting  the  largest  of  these  anchors  from  a  reasonable  gun  tube 
requires  projective  forces  (and  therefore  recoil  forces)  of  ever  1.5  million 
pounds.  There  is  no  practical  way  to  tie  the  PEA  gun  barrel  directly  to  the 
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earth  during  normal  use.  Therefore,  the  handling  of  the  large  recoil  forces 
generated  by  PEAs  is  a  significant  problem.  It  is  the  topic  of  this  report. 

It  will  be  necessary  to  develop  a  quantitative  understanding  of  the 
problem  before  looking  at  potential  solutions.  Figure  2  shows  the  general 
layout  of  an  unsupported  launch  tube,  propellant  charge,  and  anchor.  It  also 
shows  the  principal  forces  acting  on  anchor  and  launch  tube.  The  net  force 
acting  on  the  anchor  and  the  launcher  is: 


A  - 1  - 1 


F. 

■ l 


(P-Pa) 


1/  P  C 
I'l 


di 


U2 

i 


where 


subscript  i 


=  designates  anchor  when  i  =  1  and  launch 
tube  when  i  =  2 
=  net  force  (pounds) 

=  gas  pressure  in  launch  tube  (psf) 

*  ambient  pressure  (psfa) 

-  cross-sectional  area  of  launch  tube  and  piston  (sq  ft) 
=  friction  force  between  anchor  piston  and  launch 
tube  wal 1  (pounds) 

=  density  of  seawater  (slugs/cu  ft) 


A.  =  maximum  cross  sectional  area  of  anchor  or 
launch  tube  (sq  ft) 

=  drag  coefficient  (dimensionless) 

■  =  velocity  (ft/s) 


The  instantaneous  accelerations  of  the  anchor  and 


1  aunch 


tube  are: 


■A-l-2 


where 


U.  *  a.  =  F./M. 

i  i  i.  i 


U.  .  *  a.  e  acceleration 

i  i 


(ft/s2) 

*  mass  (including  virtual  miss)  (slugs) 
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The'  instantaneous  velocities  are: 
A-I-3  U. 


.  rv 

1  Jn 


where 


0 


t  =  time  (sec) 


The  distances  travelled  by  the  anchor  and  launch  tube  are  given  by: 

/t 

A- 1 -4 


X.  =  f  U.  dT 

1  Jq  1 


where 


x.  -  the  distance  traveled  by  the  anchor  or  the 
1  launch  tube  (ft) 


All  one  must  do  to  compute  the  trajectories  of  the  anchor  and  launch  tune  is 
to  integrate  equations  A-I-3  and  4  for  i  =  i  and  2  simultaneously  until 


A-I-b  xj  +  x2  =  L 

where  L  =  the  effective  length  of  the  launch  tube  (ft). 


One  must  know  the  various  coefficients  and  boundary  conditions  before  the 
equations  can  be  integrated.  In  addition,  the  pressure  must  La  known  as  a 
function  of  time  and/or  the  sum  In  the  detailed  calculations  the 

pressure  will  be  given  as  a  differential  equation  involving  pressure,  time, 
and 


for  approximate  calculations,  one  may  neglect  the  F^  and  the  fluid  dy¬ 
namic  drag.  Further,  the  pressure  may  be  assumed  constant  because  that  gives 
the  best  performance  and  because  it  is  nearly  achievable*  With  these  simpli¬ 
fications,  the  final  velocity  of  the  anchor  is 


A-I-b 


U 


1  final 


V 


2(F-THT 


"2 

k-*2 


The  first  term  in  the  radical  on  the  right-hand  side  of  Equation  A-I-b  is  the 
square  of  the  final,  anchor  velocity  that  would  result  if  the  launcher  were 
firmly  attached  to  a  very  large,  effectively  infinite,  mass.'  The  second  term 
on  the  right-hand  side  is  a  correction  term  for  the  case  where  the  launcher  is 
attached  to  a  finite  mass,,  The  product  of  these  two  thrms  is  the  final  anchor 
velocity  when  the  mass  of  the  launch  tube  is  significant  and  free  to  move. 
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The  correction  term  is  always  less  than  one.  It  increases  as  M9  increases 
relative  to  Equation  A  - 1  -6  shows  that  U.^  final  can  be  increased  by 

increasihg  P,  A,  L,  and/or  and/or  by  decreasing  . 

II.  CURRENT  PEA  DESIGNS:  To  increase  the  effective  mass  of  the  launcher 

without  increasing  its  physical  mass,  the  original  designers  added  reaction 
vessels  to  increase  the  v4rtual  mass  of  the  launcher.  Virtual  mass  is  the 
mass  of  the  water  that  gets  accelerated  ?iong  with  an  object.  The  virtual 
mass  is  added  to  the  physical  mass  of  the  object  being  accelerated,  but  only 
while  it  is  being  accelerated.  The  quantity  of  the  virtual  mass  is 
proportional  to  the  maximum  cross-sectional  area  of  the  object,  taken 
perpendicular  to  the  direction  of  acceleration.  The  proportionality  constant 
is  determined  by  the  shape  of  the  object.  Figure  3  shows  the  present  PEA 
design  and  how  the  virtual  mass  is  increased  by  adding  the  reaction  vessels. 
This  is  a  good  shape  for  increasing  the  virtual  mass  as  well  as  .increasing  the 
drag.  While  this  design  achieves  the  required  anchor  velocity,  it  is  large, 
heavy,  and  awkward.  Furthermore,  the  reaction  vessels  on  che  larger  ancho' 
launchers  occasionally  break.  The  rest  of  this  paper  is  a  description, 
analysis,  and  evaluation  <«f  several  new  concepts  for  launchers  that  do  not 
need  reaction  vessels.  ’ 

III.  PROPOSED  NEW  APPROACHES  TO  LAUNCHER  MASS  REDUCTION:  While  the 
solutions  developed  in  this  paper  also  apply  to  the  smaller  PEAs,' the  main 
emphasis  is  placed  on  the  300  KIP  PEA  because  of  its  size,  weight,  and  han-  . 
dling  problems.  The  purpose  of  these  new  approaches  is  to  accelerate  the 
anchor  to  the  required  velocity  without  making  the  launcher  large,  heavy, 
awkward,  or  prone  to  failure.  AH  of  the  new  approaches  have  one  feature  in 
common.  Instead  of  reducing  the  effect  of  the  recoil  force  by  adding  mass  to 
the  launcher,  they  all  neutralize  tlte  recoil  force  with  another  source  of. 
thrust. 

The  three  basic  options  are  shown  in  Figure  4.  The  first  option  is  the 
archetype  for  the  rest.  It  is  a  sca1ed*up  recoilless  rifle  with  the  anchor  as 
a  projectile.  Conceptually,  it  is  a  PEA  with  an  enlarged  gas  generator  and  a 
rocket  nozzle  venting  to  the  rear.  The  thru  ;t  from  the  rocket  nozzle  can  be 
tailored  to  exactly  neutralize  the  recoil  force.  This  makes  effectively 
infinite.  Since  there  is  no  recoil  displacement  to  reduce  the  efficiency  of 
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the  launcner,  the  pressure  may  be  reduced.  This,  in  turn,  allows  the  launch 
tube  to  be  thinner  and  lighter.  On  the  negative  side  is  the  need  for  a  much 
larger  and  more  expensive  gas  generator. 

The  second  option  is  derived  from  the  first  by  the  observation  that  the 
rocket  thrust  required  t,  overcome  the  recoil  force  is  exactly  equal  to  that 
required  to  propel  tne  anchor.  Therefore,  why  not  eliminate  the  launch  tube 
and  propel  the  ancnor  directly  by  the  rocket's  thrust?  The  launch  tube  is 
eliminated  completely.  There  are  other  less  obvious  advantages  to  this  op-  ■ 
tion.  The  thrust  duration  is  not  limited  by  tne  ’ength  of  a  launch  tube,  • 
because  there  is  no  launch  tube.  The  thrust  duration  is  only  limited  by  the 
standoff,  distance  between  the  sea  floor,  and  the  launch  position.  The  rocket 
thrust  can  be  reduced  because  the  required  thrust  is  inversely  proportional  to 
this  standoff  distance.  The  principal  disadvantage  is  that  the  'mass  of  the 
rocket  motor  must  be  added  to  M,,  the  mass  of  the  anchor.  This  requires  an 
even  larger  and  more  expensive  vented  gas  gene. ator  than  that  required  for  the 
recoil  less  launcher. 

The  third  option  is  similar  to  the  first  in  overall  concept.  Instead  of 
using  an  ordinary  rocket  which  vents  the  propellant  gas,  this  concept  uses  a 
water  jet.  Gas  from  a  gas  generator  forces  water  out  the  nozzle.  The  princi¬ 
pal  advantage  of  the  water  jet  is  a  smaller  (by  half),  simpler,  and  unvented 
gas  generator.  An  unvented  gas  generator  avoids  all  of  the  problems  of  the 
large  vented  gas  generators. 

3.  RECOILLESS  OPTION 

i  ‘ 

While  all  of  the  options  considered  in  this  report  are- effectively 
recoi l less, . the  term  will  be  reserved  for  the  combination  of  a  launch  tube 
with  a  rocket  which  neutralizes  the  recoil  force.  The  recoilless  option  is  a 
direct  copy  of  the  recoi  Hess  rifle  first  used  expensively  in  the  Korean  War. 
This  original  oecoilless  rifle  eliminates  the  need  for  hteavy  gun  mounts  by 
eliminating  the  recoil,  'nis  permits  light  artillery  to  he  fired  from  light 
tripod  mounts.  These  lignt  tripods  are  required  only  to  support  the  weight  of 
the  gun  and  to  facilitate  aiming;  they  do  not  have  to  absorb  any  recoil 
forces. 
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I.  RECOILLESS  RIFLF  TECHNOLOGY  APPLIED  TO  PEAs 

As  shown  in  Figure  4,  a  recoiiless  anchor  launcher  has' a  single  combus¬ 
tion  chamber  which  produces  gas  to  drive  the  anchor  via  the  piston  and  to 
power  the  rocket.  The  piston  area  and  the  nozzle  throat  area  are  preset  so 
that  the  two  forces  acting  on  the  launch  tube  are  equal  and  opposite,,  thus 
keeping  the  launcher  stationary.  As  before,  the  force  accelerating  the  anchor 


8-1  -I 


FV  *  (P  -  V  A  -  ff  -  V2  *  A!  Cdl  U* 


The  recoil  force  is: 


8-1-2 


F2  =  (R  -  V  A  -  Ff 


p  ... 
/ 


The  thrust  from  the  rocket  must  equal  F ^  if  the  launcher  is  to  be  kept 
stationary.  The  rocket  thrust,  is: 


B-I  -3 


where 


F3  =  P  At  Cf 


=  throat  area  of  ttw  nozzle  (sq  ft) 
C,  =• thrust  coefficient  (dimensionless) 


Equating  equations  B-I -2  and  3,  and  solving  for  A  gives: 


8-1-4 


A  (P-Pa)-F. 


Pfl  is  small  with  respect  to  P,  and  Ff  is  fairly  small.'  Therefore,  A^  can  be 
estimated  by: 


8-1  -h 


A  =  — 
t  C, 


The  net  force  on  the  launch  tube  for  the  recoilless  option  is: 


8-1  -b 


f2T  MP  -  Pa)  A  -  Ff  -  P  At  Cf  -  V2p  A2C0;,  Uj 
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To  do  a  detailed  calculation  of  the  performance  expected  from  any  given 
set  of  parameters,  one  must  substitute  equations  B -I - *  and  6  into  equations 
A - I -2  and  integrate  as  in  equations  A-I-3  and  4  to  the  boundary  condition  ex¬ 
pressed  in  equation  A-I-5.  To  do  this  one  needs  an  equation  (algebraic  or 
di f ferenti al )  for  P. 

If  one  assumes  that  P  is  constant  and  F g j  is  nearly  zero,  then  the 
equations  can  be  solved  analytically.  The  general  equations  for  the  anchor 
velocity  and  displacement  are: 


B-I-7 

,U1  = 

U^tanhU/i)  and  • 

R-I-8 

X1 

ul^c  loge(cosh(t/T)). 

where 

ui 

1*> 

terminal  velocity  of  the  anchor  if  the  driving  force 

were  maintained  indefinitely.  When  Uj  equ  Uj«» 

then  F^  in  equation  B-I-l  =0. 

B-I-9 

u, 

1“> 

(P-Pa)  A  -  Ff 

1/2pAlCdl  ■ 

B-I-1U 

T  = 

a  characteristic  time 

T  — 

M1  Ul» 

ip-VA-Ff 

The  time,  t,  can  be  eliminated  from  equations  B-I-7  and  8  to  give  as  a 
function  of  x. 


8-1-11  ,'Uj  =  U1(o  1  -  exp  (-2x/Uia>r) 

where  exp  (3)  =  e^ 

The  final  velocity  can  be  calculated  by  substituting  L  for  x  in  equation 
B-l-1 1 . 

These  equations  are  useful  for  evaluating  proposed  systems  and  for  com¬ 
paring  different,  types  of  systems,  but  they  don’t. describe  the  gas  generator 
necessary  to  provide  the  constant  thrust.  No  real  gas  generator  provides  a 
perfectly  constant  pressure;  therefore,  more  detailed  calculations  are  needed 
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to  give  accurate  results.  These  detailed  calculations  will  give  answers  for  a 
real  gas  generator,  but  more  important,  they  will  be  a  powerful  design  tool 
for  the  development  of  the  gas  generator. 

II.  INTERIOR  BALLISTICS:  This  technology  area  gives  the  means  to  calc¬ 
ulate  the  changing  interior  pressure,  P,  and  the  means  to  design  an  effective 
propellant  system.  The  first  equation  needed  is  the  equation  of  state,  winch 
relates  the  pressure,  volume,  mass,  temperature,  and  molecular  weight  of  the 
gas.  The  ideal  gas  equation  is  usually  used  for  more  moderate  pressure,  but 
the  high  pressures  used  in  the  PEAs  need  an  equation  of  state  such  as  the 
modified  van  der  Waals  equation  of  state. 


B-ii-i  P  '('' -  »«g) 

b  =  covolume,  a  measure  of  the  finite  volume  of  the  gas 
molecules 

MW  =  molecular  weight  of  the  gas  (slugs/slug-mole) 

V  =  volume  available  to  the  gas  (cu  ft) 

■  o 

R  =  universal  gas  constant  (49,709  ft-lb/slug  mole-  R) 

o 

T  =  absolute  temperature  of  the  gas  (  R) 
m  =  mass  of  gas  (slugs) 

MW  and  b  are  parameters  whose  values  depend  on  the  composition  of  the  gas;  R 
is  a  universal  constant;  P,  V,  m,  and  T  are  variables.  To  calculate  instant-* 
aneous  values  of  P,  it  is  necessary  to  know  the  two  parameters,  MW  and  b,  the 
universal  gas  constant ,  R ,  and  the  variables  V,  m,  and  T.  The  remainder  of 
this  section  covers  the  calculations  of  these  variables. 

The  process  starts  with  the  burning  of  the  propellant  to  produce  the 
gas.  The  surface  of  the  propellant  regresses  everywhere  perpendicular  to  the 
burning  surface,  figure  5  shows  a  propellant,  grain  with  a  single  perforation 
(monoperf),  and  Figure  6  shows  a  grain  with  seven  perforations  (seven  perf). 
Figures  7  and  B  show  a  rocket  propellant  grain  of  a  design  that  might  be 
appropriate  to  a  recoil  less  or  direct  rocket  PEA.  The  drawings  show  how  the 
grains  regress  as  they  burn;  the  grain  is  shown  unburned  and  partially 
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burned.  The  following  equations  give  the  surfaces  of  the  three  types  of 
grains  as  functions  of  their  original  dimensions  and  the  burned  distance,  z. 
These  functions  are  needed  to  calculate  the  rate  of  gas  production. 

B-II-2  S  =  *  (OD  +  ID  )  {  V?  (OD  -  ID  }  +  H  -  4z} 

m  .  •  o  o  1  c  o  o  o  J 

S 

where  m  •  . =  burning  surface  of  one  monoperf  propellant 

grain  (sq.  ft.) 

0Dq  =  original  outside  diameter  of  grains  (ft.) 

Ig  =  original  inside  diameter  of  perforation  (ft) 

H  =  original  length  of  grain  (ft) 

2°  =  distance  burned,  perpendicular  to  grain's 

surface  (ft) 

This  monoperf  grain  is  totally  consumed  whenever  z  =  Hq/2  or  z  =  (00  q  -  1Dq)/4, 
whichever  comes  first.  The  equation  for  the  seven  perf  grain  is: 

B-II-3  S  =  n  {  V?  (OD  -7 ID  )+H  (OD  +71D)  +2z(6H  -20D  -14ID  )-3bz2} 

S  lt0,000  0  0  0  0 

where  S?  =  burning  surface  of  one  seven  perf  propellant 

grain  (sq  ft) 

OD  -  '31 D 

This  equation  for  5  ,  is  valid  while  z  >  — 2 — ^ — 2.  .  when  z  reaches  this 

limit,  7/8  of  the  grain  is  burned.  The  remaining  part  of  the  grain  burns  in 
two  modes.  It  takes  about  six  equations  to  calculate  the  surface  of  these 
afterburns,  these  equations  are  not  shown  here,  but  they  are  in  the  computer 
program  listed  in  Appendix  A.  The  equation  for  the  rocket  propellant  grain 
is: 

B-II-4  S  '  ■  (H  2w  -  4z)  116  h  *  8  ID  f  +  ir  (ID  +8w  +  2z)} 

r  o  o  o  o 

-2t»(3w  -  z)(w  ♦  2)  *  41  Dq { <i> ( *  4wV  -  2w  cos 
Sr  =  instantaneous  surface  of  a  rocket  type  grain  (sq  ft) 
hQ  *  length  of  straight  sides  of  the  «fokes  (feet) 
w  =  web  of  propellant  grain  (feet). 

$  =  arcsine  (2w/IDq)  (radians) 
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The  next  item  needed  is  the  propellant's  regression  rate,  z  .  ■ 

B-II-5  z  =  'B  Pn  , 

where  B  =  burning  rate  coefficient  (ft/s  psf  n) 

n  =  burning  rate  exponent  (dimensionless) 

Both  B  and  n  are  determined  by  the  composition  of  the  propellant.  The  overall 
burning  rate  mb  is  given  by: 

•  i  « 

m,  =  N  S  p  z 
b  p 

N,  =  number  of  propellant  grains  (dimensionless) 

S  =  Sm ,■  Ss,  or  Sr  as  appropriate  (sq  ft) 
p  =  density  of  the  propellant  (slugs/cu  ft) 

p 

In  both  the  recoil  less  and  the  rocket  launched  PEAs,  most  of  the  propellant 
gas  goes  out  the  nozzle.  The  equation  for  the  exhaust  mass  flow  rat,e  is: 

B-II-7  m  =  P  A  C 

e  t  e  ■ 

where  At  -  area  of  the  throat  of  the  nozzle  (sq  ft) 

Cg  =  discharge  coefficient  of  propellant  gas  (slugs/pound) 
Typical  values  for  the  discharge  coefficient  are  .0002  slugs/pbund-seconds. 


The  net 

rate  of  accumulation 

of  propellant  gas  is: 

B-H-8 

• 

m  * 

N  S  Pp  B  Pn  -  P  Af  Ce 

where 

m(0 )  * 

P(0)V(0)  MW 

b  P(0)  +  R.T(O) 

The  next  item  to  be  calculated  is  the  volume,  V,  available  to  the  gas.  This 
volume  comes  from  three  sources;  the  initial  or  free  volume,  Vf;  the  rate  of 
increase  of  volume  due  to  motion  of  the  piston  which  is  equal  to  (U^  ♦  Ug)  A  ; 
and  the  burning  of  the  propellant  grains.  The  gas  volume,  V,  is  given  as  a 
derivative  with  respect  to  time.  The  boundary  condition  at  time  *  0  is  . 


B-II-6 

where 
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B-II-9  V  =  (U,  +  U2)A  +  N  S  B  Pn 

and  V(0)  =  V^.  =  initial  volume  available  to  gas  (cu  ft) 

The  next  state  function  to  be  calculated  is  the  temperature,  T.  The 
changes  in  temperature  will  be  rather  small  in  all  of  the  configurations  being 
considered  here,  because  the  propellants  burn  to  produce  gas  at  constant  temp¬ 
erature.  This  constant  flame  temperature  and  the  high  burning  rates  tend  _o 
maintain  the  temperature  nearly  constant.  The  combination  of  a  high  through¬ 
put  and  a  short  action  time  means  that  heat  transfer  is  minimal.  The  only 
significant  heat  loss  comes  from  the  expansion  work  done  by  the  gas.  The  time 
derivative  of  the  temperature  is: 

B-II-10  T  =  (T,  -  T)  m. /m  -  P  V  /(c  J  m) 

f  b  .  u  p 

cp  =  specific  heat  (Btu/slug-°R ) 

J  =  mechanical  equivalent  of  heat  (778.3  ft-lb/Btu) 

r  ]  -  _!_  JL 
p  Y-l  MW  •  • 

Y  =  ratio  of  specific  heats  (dimensionless) 


! 

T  (0) 

,-Y 

, 

■ 

* ' 

V 

u 

■  tu(*  U2)A 

* 

i  ? 

Substitute  Equations 

1 

B-l.-l  and  6  into  Equation 

A-I-2  gives: 

8-II-11 

• 

U 

1  *  {  (P  -  Pfl)  A  -  Ff  -  l/2  p 

Ai  cdi  ui  i 7ni 

» 

B - II  - 1 2 

•  i 

u 

2  '  *  %)  A  -  Ff  -  1/2  p 

A2  Cd2  ^2  *  P  At  Cf^  /M2 

*• 

• 

Taking  the 

deriviti v 

e  of  equation  A  - 1 -4  gives: 

* 

i 

B-ll-13 

• 

X 

1  =  U1 

- 

B*tl-14 

* 

X 

2  *  U2 

a 

,11  . 

1 

.  _  ■ 

•  v  •. 

.  **  *.  v  ♦. .  •,  •:  v  % 

...  .  ... 

■  ' 

.  ' 

'  1  **•*_. 
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All  of  the  equations  for  the  recoilless  launcher  are  finally  assembler!. 
The  first  step  in  calculating  the  performance  of  a  specific  recoilless  design 
is  to  determine  all  of  the  parameters  and  boundary  conditions.  The  next  step 
is  to  select  an  appropriate  burning  surface  equation;  equation  B-II-2,  3,  or 
4,  or  some  other  equation  if  a  different  grain  is  used.  The  burning  surface 
and  equation  of  state,  equation  B-II-1,  are,  auxiliary  equations.  Equations 
B-II-8  through  B-II-14  are  all  first  order  differential  equations  that  must  be 
integrated  simultaneously  until  the  final  boundary  condition  is  met, 

(xj  +.x^  =  L).  Because  of  the  mixed  nature  of  these  equations,  it  will  be 
necessary  to  use  numerical  techniques  to  integrate  them.  If  the  reader  is 
planning  to  write  a  computer  program  to  carry  out  these  integrations,  he 
should  be  aware  of  several  conventions  used  in  the  derivation  of  the  equa¬ 
tions:  The  sign  of  is  positive  when  it  is  going  forward,  the  desired 

direction.  Similarly,  the  sign  of  U,,  is  considered  positive  when  it  is  going 
rearward,  its  natural  direction  of  motion.'  might  become  negative  if  the 
rocket  thrust  is  higher  than  the  recoil  force.  It  is  extremely  unlikely 
that  will  ever  become  negative,  but  the  possibility  that  U2  can  become 

negative  may  cause  a  problem  in  equation  8- 11-12.  All  of  the  terms  in  this 
equation  have  intrinsic  directions  indicated  by  their  signs.  Tbe  fluid  drag  . 
term  {  V2  pC^^A  .1) .  )is  always  opposite  in  sign  to  the  velocity  The 
equations  as  given  are  correct  so  long  as  the  velocity,  U^,  is  positive.  This 
can  be  handled  in  a  computer  program  by  multipying  this  term  by  a  function 
SIGNUM(Ug).  This  function  has  a  value  of  *1  or  -1  as  as  tl^  is  positive  or 
negative. 

Appendix,  A  is  a  copy  of  a  computer  program  like  the  one  just  described, 
and  a  rub  showing  a  typical  set  of  results.  This  program  was  encoded  several 
months  before  the  detai led' analysi s  shown  in  this  report.  This  explains  the 
difference  between  them.  For  example,  the  program  yses  the  ideal  gas  equation 
of  state  , instead  of  the  modified  van  der  Waals  equation  of  state.  Also  the 
program  does  not  have  a  term  for  fluid  dynamic  drag  on  the  launcher,'  because 
the  launcher  velocity  is  intended  to  be  so  low  that  the  drag  is  insignificant. 

The  calculations  show  that  the  anchor  can  be  accelerated  to  300  ft/s  by  a 
recoil  less  device  operating  at  20,000  psi.  The  propellant  charge  assumed  in 
this  calculation  consisted  of  115  seven  perf  grains  (Figure  6).  This  is, not' 
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the  usual  way  to  fuel  a  solid  propellant  rocket.  Most  solid  propellant  rocket 
motors  have  one  or  two  grains  that  are  nearly  as  large  as  the  motor.  The 
grains  are  also  of  such  a  shape  that  they  stay  in  place  with  a  very  high' 
reliability.  This  is  not  the  case  with  relatively  small  grains,  they  can  be 
blown  through  the  large  nozzles  required  by  the  PEAs.  It  is  not  clear  that 
the  techniques  for  grain  retention  that  work  in  the  much  smaller  recoiliess 
rifles  will  work  in  the  PEA  combustion  chambers.  If  this  potential  problem 
becomes  serious,  then  it  will  be  necessary  to  use  rocket-like  propellant 
grains  as  described  in  Figure  7,  if  they  can  be  made  to  sustain  the  high  acce¬ 
leration  load.  It  seems  clear  that  the  development  of  the  propellant  assembly 
would  be  an  expensive  and  chancy  undertaking. 

C.  DIRECT  ROCKET 

The  recoil  less  PEA  launcher  system  described  in’ the  previous  section 
appears  to  be  a  theoretically  workable  solution  to  the  problems  involved  in 
launching  large  PEAs,  if  the  propellant  problems  can  be  solved.  Also,  it  is 
rather  sophisticated.  This  level  of  sophistication  is  accepted  in  light 
artillery  because  it  can  give  high  velocity  and  accuracy  without  recoil  in  a 
■hand  held  weapon'.  PEAs  are  not  hand  held  and  they  do  not  require  high  accur¬ 
acy.  High  velocity  and  adequate  accuracy  can' be  acnieved  with  unguided  direct 
rocket  propulsion,  but  the  direct  rocket  shares  the  recoilless  PEA’s  propel¬ 
lant  problems.  With  the  direct  rocket,  the  gun  tube  and  piston  elimin- 
ated.  '  ,  •  , 

I.  DESCRIPTION  OF  DIRECT  ROCKET  PROPELLED  PEA: 

A  direct  rocket  powered  PEA  i.s  shown  in  Figure  4.  The  rocket  motor  is 
attached  directly  to  the  anchor.  The  thrust  from  the  rocket  motor,  less  the 
fluid  dynamic  drag  on  the  assembly,  accelerates  the  anchor  and  the  rocket 
motor.  The  analysis  of  the  forces  is  simpler  than  for  the  original  PEA  system 
and  the  recoil  less  system,  because  the  whole  apparatus  moves  as  a  unit. 

The  "up  front"  advantages  of  the  direct  rocket  propelled  PEAs  are  the 
reduced  weight  and  complexity  resulting  from  the  elimination  of  the  launch 
tube  and  piston.  There  is  another,  less  obvious,  advantage  to  the  direct 
rocket  PEAs.  The  duration  of  the  thrust  is  not  tied  to  the  barrel  length, 
because  there  is  no  barrel.  The  rocket  can  fire  for  a  somewhat  longer  time  at 
a  lower  thrust  level.  If  desirable,  the  thrust  can  be  continued  through 
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penetration  of  the  sea  floor.  The  direct  rocket  PEAs  also  have  several  disad¬ 
vantages  relative  to  the  original  and  the  recoil  less  PEAs.  In  the  original 
and  recoilless  systems,  only  the  anchor  and  the  piston  are  accelerated  for¬ 
ward.  In  the  direct  rocket  PEA,  everthing  is  accelerated  forward.  This 
substantially  increases  the  drag  and  the  propelled  mass,  thus  requiring  an 
increased  thrust.  Another  disadvantage  of  the  direct  rocket  is  that  the 
propellant  grain(s)  is  exposed  to  the  same  high  acceleration  and  gas  flow  as 
the  anchor  (250  g's).  This  could  cause  the  propellant  to  break  up  and  the 
motor  to  fail  or  rupture. 

II.  ANALYSIS  OF  DIRECT  ROCKET  PEA: 

This  analysis  is  less  wordy  because  it  closely  parallels  the  analysis  of 
the  recoilless  PEA.  It  starts  with  the  forces  on  the  body.  The  subscript,  3, 
will  be  used  to  distinguish  variables  relating  to  the  di rect  rocket  PEA.  The 
net  force  on  the  body  is: 

c-n-l  F3  ■  P  \  Cf  -  »  u3  ft3  cd3 

All  of  the  terms  have  been  defined  earlier.  The  terminal  velocity,  U^,  is 
the  velocity  at  which  F3  =  \  It  can  be  calculated  by  setting  F-j  equal  to 
zero  in  Equation  C-II-1,  and  solving  for  the  terminal  velocity. 

C-1I-2  = 

The  rate  of  change  of  the  mass  of  the  vehicle  is.: 

C-tI-3  H-  *  -  P  A  C 

3  ■  t  e 

The  initial  mass  of  the  vehicle,  Including  propellant,  is  M3(Q).  The  velocity 
of  the  vehicle  is  given  by  the  differential  equation: 

C-TI-4  '  .  U3  *  F3/M3 
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The  position  of  the  vehicle  is  given  by  the  differential  equation: 

C-II-5  X3  =1)3 

The  initial  values  of  U3  and  X3  are  zero.  All  that  is  necessary  to  integrate 
the  equation  is  a  way  to  calculate  the  instantaneous  value  of  P.  This  can  be 
done  with  most  of  the  same  equations  used  for  the  recoilless  PEA.  Equations 
B-lI-1  to  8-1 1 -8  can  be  used  as  they  are.  The  rate  of  change  of  the  gas  vol¬ 
ume  is  different  because  there  is  no  piston. 

C-II-6  V3  =  N  S  8  Pn 

Equation  B-II-10  is  not  needed  because  the  temperature  remains  constant. 

As  with  the  recoil  less  PEA  calculations,  the  first  step  is' to  get  values 
for  the  boundary  conditions  and  parameters.  A  burning  surface  equation  must 
be  selected  from  B-II-2,  3,  or  4;  or  some  combinations  of  them;or  some  other 
equation  as  needed.  The  equation  of  state,  equation  B-II-1,  is  needed,  and 
the  differential  equations  B-II-8,  C-II-4,  5,  and  6  are  needed.  Again,  as 
with  the  recoil  less  PEA,  these  equations  must  be  solved  by  numerical  techni¬ 
ques.  An  example  program  is  shown  in  Appendix  B  along  with  a  sample  calcu¬ 
lation. 

The  differential  equations  can  be  solved  in  closed  form  if  the  pressure 
is  constant.'  The  velocity  as  a  function  of  time  is: 


C-II-7 


where 


U3  =  U ^  tanh 


>«,c.  1 


C  =  C-/C  *  velocity  of  rocket  exhaust  (ft/s) 
f  e  • 


D.  MATER  JET  PEAs 

While  explaining  the  details  of  the  assignment  to  the  author,  the  spon* 
sors  mentioned  an  unusual  propulsion  scheme  proposed  at  Port  Hueneme.  That 
.proposal  was  for  a  gas  driven  water  jet  to  neutralize  the  recoil  force.  Our 
understanding  was  that  the  gas  was  to  be  supplied  from  compressed  gas  (air) 
cylinders.  To  get  the  needed  re$ults  would  require  nearly  SO ,000  standard 
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cubic  feet  of  gas  compressed  to  40.000  psi .  Half  of  this  compressed  gas  would 
have  to  be  transferred  to  the  water  chamber  in  40  milliseconds.  This  would 
require  a  very  impressive  valve.  While  it  seems  very  unlikely  that  the  job 
can  be  done  with  compressed  gas,  the  picture  is  quite  different  if  the  driving 
gas  is  generated  by  burning  a  charge  of  gun  propellant.  Instead  of  nearly  twt 
tons  of  very  high  pressure  air,  not  to  mention  the  tankage,  the  job  can  be 
done  with  only  200  pounds  of  gun  propellant. 

I.  WATER  JET  PEA  DESCRIPTION 

Figure  4  gives  a  conceptual  drawing  of  a  water  jet  PEA.  A  water  jet  PEA, 
operating  at  20,000  psi V  (2,880,000  psf),  will  have  a  six  foot  long. launch 
tube  with  an  inside  diameter  of  ten  inches  (.8333  ft).  This  would  be  connect¬ 
ed  to  the  end  of  water  tube  that  is  5.5  feet  long  and  a  2.12  foot  inside  diam¬ 
eter.  Its  shape  will  be  .cylindrical  with  hemispherical  ends.  The  free  end  of 
the  5.5  foot  diameter  tube  is  terminated  with  a  Convergent  nozzle  with  an  exit 
area  of  .25  square  feet  (ID  =  .564  ft).  The  propellant  charge  and  the  igniter  ' 
are  stored  in  the  end  of  the  water  tube  near  the  launch  tube;  the  drive  piston 
fits  into  the  .833  foot  diameter  tube;  and  the  seawater  goes  into  the  2.12 
foot  inside  diameter  water  tube.  The  water  tube  is  dry  during  storage  and 
handling.  It  fills  with  seawater  while  being  lowered  into  the  sea. 

The  propellant  charge  is  ignited  when  the  assembly  is  lowered  to  its 
deployment  depth.  The  burning  propellant  causes  the  pressure  to  rise  quickly 
to  a  constant  level  of  about  20,000  psi  internal.  The  high  pressure  causes 
.the  water  to  accelerate  and  stream  out  of  the  nozzle.  After  a  few  millisec¬ 
onds,  the  exit  velocity  stabilizes  at  about  1,700  feet  per  second.  Near  the 
end  of  the  water  expulsion,  there  is  another  acceleration  as  the  nozzle 
empties.  Figure  10  is  a  graph  of  the  water  velocity  as  afunction  of-  time. 

The  higher  velocity  trace  represents  the  velocity  of  the  water  leaving  the 
nozzle.  The  lower  velocity  trace  represents  the  velocity  of  the  gas-water 
interface. 

II.  ANALYSIS  OF  THE  WATER  JET  PEA- 

The  equations  needed  to  calculate  the  performance  of  the  water  jet  PEA 
are  the  same  as  those  used  for  the  recoil  less  PEA,  with  one  exception.  That 


1  A  later  section  will  cover  water  jet  PSAs  operating  at  Other  pressures. 

•  .  \  '* 

'  \ 
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exception  is  equation  B-I-3,  which  gives  the  thrust  from  the  rocket.  The 
thrust  from  the  water  jet  is: 

D-n-i 

where  U^,=  exit  velocity  of  water  (ft/s) 

U4  =  velocity  of  water  at  interface  with  drive  gas  (ft/s) 

L  =  length  of  the  water  tube  (feet) 

The  definition  of  the  terms  is  given  p i ctori ally  in  figure  9.  The  detailed 
calculation  is  somewhat  lengthy  and  complex,  and  it  will  not  be  derived 
here.  The  derivation,  along  with  a  computer  model  and  sample  problem,  is 
given  in  Appendix  C.  If  one  assumes  a  constant  gas  pressure  of  2Q,.QIQG  psi , 
the  following  events  occur  in  sequence.  Figure  10  should  be  referenced  during 
the  explanation.  Initially,  the  water  is  at  rest,  but  it  starts  to  accelerate 
as  the  pressure  is  applied.  The  initial  acceleration  of  the  main  body  of 
water  is  about  240,000  ft/s2  but  it  drops  to  zero  in  about  four  millisec¬ 
onds.  The  average  acceleration  over  the  four  milliseconds  is  42,500  ft/s2. 

The  initial  high  acceleration  is  caused  by  the  nigh  pressure  gradient  across 
the  constant  area  length  of  the  water  column.  As  the  velocity  increases,  the 
pressure  gradient  across  the  constant  area  length  decreases  until  the  accel¬ 
eration  reaches  zero.  At  that  time,  all  of  the  pressure  gradient  is  across 
the  nozzle  where  the  speed  increases  by  a  factor  of  ten.  A  particle  of  water 
in  the  constant  area  section  is  traveling  at  about  170  ft/s,  so  long  as  it 
remains  in  the  constant  area  section.  When  it  reaches  the  nozzle  it  accel¬ 
erates  in  one' mi  1 1 i second  to  1,700  ft/s,  an  acceleration  of  1,530,000  ft/s^. 

Without  detail,  the  operating  sequence  is:  Four  milliseconds  acceler¬ 
ating  the  main  body  of  water  to  its  steady  state  speed  of  170  ft/s;  35  milli¬ 
seconds  at  the  steady  speed  of  170  ft/s;  and  one  millisecond  of  very  high 
acceleration  as  the  nozzle  empties.  During  this  whole  process,  the  water 
leaving  the  nozzle  is  traveling  ten  times  faster  than  the  water  in  the  main 
body.  During  the  main  steady  portion,  the  exit  velocity  is  1,700  ft/s.  The 
exit  velocity  can  easily  be  calculated  fur  the  steady  flow  by  equating  the 
total  heads  at  the  inlet  and  the  exit  of  the  nozzle  along  with  the  continuity 
equation.  • 
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D-II-2  ■  P  +V2  p  Uj  =  +1/2  p  U| 

whe^e  =  velocity  of  the  gas -water  interface  (ft/s) 

U5  =  water  velocity  at  exit  from  nozzle  (ft/s) 

The  continuity  equation  is: 

P  Awt  U4  =  p  At  U.5 

A  =  cross  sectional  area  of  straight  portion  of  the 
w  water  tube  (sq  ft) 


D-II-3 

where 


D-II-5 


„  2  \  (MV 

5  "  1 


If  the  forces  are  balanced  so  that  the  launcher  remains  stationary,  then 
the  motion  of  the  anchor  will  be  described  by  equations  B-I-7  to  10.  The 
following  equation  must  be  satisfied  to  have  the  forces  balanceJ. 


0-11-6 


A  *  A 
t  wt 


(P-Pa')  A  *  Ff 


V^Ff 


Since  Is  small  compared  to  (P-P^)  A,  a  good  approximation  for  .  is 


0-II-7 


Assuming  the  same  2,88U,000  psf  (20,000  psi)  that  was  used  for  the 
recrilless  and  direct  rocket  examples,  the  water  jet  exit  velocity  is  !,7U4 
ft/*.  With'  a  nozzle  throat  area  A^,.  of  .25  sq  ft,  the  thrust  would  be  1.5  x  10* 
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pounds.  The  mass  flow  out  the  nozzle  will  be  852  slugs/s  (27,40U  pounds/s) 
for  .04  seconds.  The  total  mass  of  water  is  34  slugs  (weighing  1,100 
pounds).  To  hold  this  much  water,  the  internal  volume  of  the  water  tube  nust 
be  17  cu.  ft..  For  example,  the  tube  could  have  an  inside  diameter  of  2.12 
feet  and  a  length  of  b.b.foet.  The  launcher  also  requires  the  piston  tube 
which  is  six  feet  long  and  ten  inches  in  inside  diameter. 

E.  COMPARISON  OF  THE  APPROACHES 

After  having  read  this  enthusiastic  description  of  these  several  ways  of 
reliably  launching  large  PEAs,  it  is  time  to  examine  and  compare  than  critic¬ 
ally.  They  will  be  compared  on  the  basis  of  weight,  size,  and  probability  of 
successful  development.  Finally,  a  recommendation  will  be  made  as  to  which 
option  should  be  developed. 

I.  SIZE  ANU  WEIGHT  COMPARISON: 

The  size  and  weight  comparisons  are  taken  together  because  they  are  based 
on  the  same  parameters.  The  launch  tube  comparison  will  be  done  first,  be¬ 
cause  it.  is  common  to  the  recoilless  and  the  water  jet  options.'  The  key 

relationship  for  the  launch,  tubes  is  the  requirement  that  they  deliver  a  force 

* 

of  1.44  x  10  pounds  over  a  stroke  of  six  feet.  Therefore,  the  inside  diam¬ 
eter  of  the  launch  tube  is 


Note  that  inside -diameter  of  the  launch  tube  is  given  in  feet  and  the  pressure 
is  given  in  psf.  Figure  II  is  a  graph  of  the  launch  tube  inside  diarneter  as  a 
function  of  pressure.  In  order  to  determine  the  wall  thickness,' one  needs  the' 
yield. strength  of  the  \  terial  of  construction,  Because  the  pressures  will  be 
high,  it  will  be  necessary  to  use  thick  wall  equations  for  cylinders  to  calcu¬ 
late  the  wall  thickness.  This  equation  Is: 

£-1-2  P  *  S  —  or  .1  * 

where  $  *  yield  strength  of  material  (same  units  as  pressure) 

y  '  , 

and  a  00.  /  ID. 

'  .  D  D  •  . 
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Since 


00u  =  ID  +  2t 
0  b 


where  t  =  wall  thickness. 

Equation  E-I-2  can  be  rearranged  to 

VS  -* 

■  y 


E-I-3 


ID, 


t  = 


The  weight  of  the  steel  tube  launch  is  given  by 


E-I-4 


Ob  Lb  I  '002| 


I1)2b> 


where 


=  length  of  the  launch  tube  cylinder 
cb  =  density  of  structural  material  of  launch  tube. 
DC',  =  outside  diameter  of  launch  tube  {ft ) 

I; 

-  inside  diameter,  of  launch  tube  (ft) 


3ne  rust  be  careful  to  match  all  of  the  units  in  these  equations.  If  Equation 
E-I-l  and  2  are  substituted  into  Equation  E-I-4  and  simplified;  one  gets 


E-I-5 


2  F. 

=  Pfj  ~5  7YT 

y 


Figure  12  is  a  plot  of  as  a  function  of  5^  at  several  values  of  P  with  L 
fixed  at  six  feet  (72  inches).  Examining  Figure  12  shows  that  the  weight 
decreases  as  the  yield  stress  increases  and  as  the  pressure  decreases.  If  one 
reduces  the  weight  of  the  launch  tube  by  decreasing  the  pressure,  one  must 
increase  the  inside  diameter  according  to  Equation  E-I-l,  At  this  time  it  is 
not  clear  what,  if  any,  negative  effects  are  related  to  large  launch  tube, 
diameters. 
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The  launch  tube  is  the  smaller  of  the  two  main  components  of  the  proposed 
PEAs,  and  the  calculation  of  its,  size  and  weight  is  straightforward.  The  size 
and  weight  of  the  water  tube  for  the  water  jet  PEA  can  also  be  calculated  by 
strai ghtforward  analysis.  Unfortunately,  the  rocket  motors  of  the  recoil  less 
and  the  direct  rockets  are  not  so  simply  handled.  They  require  huge  thrusts 
ranging  from  1.5  to  3.0  million  pounds  and  an  extremely  short  burn  time  of 
.040  seconds.  The  volumes  of  main  stream  service  rockets  are  proportional  to 
their  total  impulse,  1^.  .  The  approximate  volume  of  a  main  stream  rocket  can 
be  calculated  by: 

E-I-6  V  =  6.2  x  10‘5  i' 

where  V  =  volume  of  a  mainstream  rocket  (cu  ft) 

Since  the  total  impul se  ;requi red  for  the  PEA  rocket  motor's  range  from  58,000 
to  80,000  pound-seconds,  the  estimated  volume  ranges  from  3.6  to  5.0  cubic 
feet.  However,  it  would  be  totally  impossible  to  make  a  rocket  of  this  volume 
deliver  1.5  to  3.0  million  pounds  of  thrust.  To  provide  the  very  large  thrust 
would  require  two  or  three  hundred  thousand  square  inches  of  burning  surface. 
The  propellant  would  have  to  be  divided  into  many  small  pieces,  and  they  must 
be  supported  in  a  dispersed  array  to  allow  room  for  6,800  pounds  per  second  of 
gas  to  flow  through  the  array  to  the  nozzle.  This  would  cause  a  large 
increase  in  the  volume  of  the  PEA  rockets  over  that  predicted  from  their  total 
impulse.  It  is  difficult  to  predict  how  large  the  volume  increase'  would  have 
to  be,  because  the  PEA  rockets  would  be  far  out  of  the  main  stream  of  r  ;ket  ■ 
motors. 

The  water  jet  must  produce  a  thrust,  F5,  of  1.44  x  10^  pounds  of  thrust 
for  .040  seconds.  Starting  with  Bernoulli's  Equation,  Equation  D-II-2,  and 
the  thrust  Equation  D-I 1-5,  one  can  derive  an  equation  for  the  mass  flow  of 
water  coming  out  from  the  nozzle.  This  equation  follows:  - 

E-I-7  -  -  -  *  •_  c  \j - - -  see  Equation  £-1-8 

w  ‘  r  5  T  P  -  P  -  Fr /A  . 

a  5'  wt 

Since  , *  ,  „  ,  r  .  , 

m  x  t.  =  m  .and  F,  t  =  l.., 
w  w*  ,  3  t3 
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We  can  multiply  both  sides  by  t,  and  calculate  the  weight  of  the  water. 


E-T-8 


Dividing  both  sides  by  pw  gives  the  volume  of  the  water  that  will  be  ejected 
during  the  launching. 


•5t 


E-i-9 


V  = 


{ 


2  p  (P  -  P 
w  .  •  a 


W 


The  dimensions  of  the  water  tube  can  now  be  calculated.  First,  the 
inside  diameter,  I03  ,  must  be  selected.  The  water  tubes  will  have  cylinders 
with  hemispherical  caps.  Given  this  shape  and  the  internal  volume,  one  can 
calculate  the  length,  L-j  . 

4  V 

.E-I-10  L  *  ' — —  *  lD,/3 

3  ,102 


Figure  13  is  a  plot  of  water  tube  length  as  a  function  of  10 3  agd  pressure. 
Because  the  walls  are  moderately  thick,  it  will  be  necessary  to  use  the  thick 
wall  equations.  The  thick  wall  equation  for  cylinders  was  given  earlier. 
Equation  E-I-2  and  3.  The  equation  for  spheres  is  very  similar 


E- 1 -11 


P 


a*  -  1 

l.b 


Since  u  is  the  ratio  of  the  outside  diameter  to  the  inside  diameter,  the 
outside  diameter  of  the  spheres  is 


E-I-12 
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The  weight  of  the  water  tube  is 


E-I-13 


wt 


=  P 


wt 


[* 


(L, 


-  ID3) 


(002 


cyl 


-  ID23 


+  £  (0Q3 

6  v  sph 


-  I03> 


Using  Equations  E-I-2,  12,  and  13,  one  can  calculate  the  weight  of  the  water 
tube.  Besides  the  equations,  cr.e  needs  the  weight  density  of  the  material, 
the  length  and  inside  diameter,  the  pressure,  and  the  yield  stress  of  the 
material . 


E-I-14 


wt 


pwt  T  P  ID3 


2(L  -  ID  )  ID 

S  -  2  P  +  S  -  l.bP 

.  y  y 


The  most  likely  material  for  the  water  tube  and  launch  tube  is  steel.  It  has 
a  weight  density  of  492  pounds  per  cubic  foot.  Composite  materials,  such  as 
fibre  glass-epoxy  would  permit  lower  weights. 


Appendix  D  is  a  copy  of  a  BASIC  computer  program  that  computes  all  of  the 
dimensions  and  weights  of  the  launch  tube  and  water  tube.  The  program  com¬ 
putes  all  of  these  values  for  a  matrix  of  pressure,  inside  diameter,  and  yield 
stress  of  material  of  construction.  The  limits  and  increments  of  the  matrix 
are  specified  by  the  user.  Besides  the  program.  Appendix  D  includes  several 
output  sheets  each  of  which  lists  the  results  for  one  inside  diameter,  a  range 
of  pressure,  and  a  range  of  yield  stresses.  Figure  11  is  a  graph  of  the 
launch  tube  inside  diameter  as  a  function  of  pressure;  Figure  12  is  graph  of 
of  the  launch  tube  weigth  as  a  function  of  yield  stress  and  pressure.  Figure 
13  is  a  graph  of  the  length  of  the  water  tube  as  a  function  of  its  inside 
diameter  and  the  working  pressure.  Figures  14  to  26  are  plots  showing  the 
weight  of  the  whole  assembly  (launch  tube  and  water  tube).  Each  figure  has  a 
family  of  curves,  one  for  each  value  of  the  yield  stress.  One  set  of  graphs. 
Figures  14  to  21  plot,  the  system  weight  versus  the  pressure.  There  is  one 
graph  for  each  of  the  eight  inside  diameters  considered.  Figures  22  to  26 
graph  :he  system  weight  versus  the  inside  diameter.  There  are  five  graphs, 
one  for  each  pressure  considered.  With  the  aid  of  these  graphs,  one  can 
select  a  configuration  which  most  closely  meets  the  needs  of  the  users..  One 
must  keep  in  mind  that  al though  , these  analyses  are  very  thorough,  they  are 
"first  cuts."  One  shortcoming  is  the  fact  that  the  sizes  given  do  not  include 
space  for  the  propellant.  , 
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F.  CONCLUSIONS  AND  RECOMMENDATIONS 

•  Upon  first  examination,  all  three  of  the  PEA  launchers  considered  in  this 
paper  look  workable.  However,  an  examination  of  the  requirements  for  the 
rocket  motors  needed  for  the  recoil  less  and  the  direct  rocket  PEAs,  shows  a 
huge  rocket  engineering  job.  The  recoilless  PEA  requires  a  total  impulse  of 
57,600  pound-seconds.  A  typical  service  rocket  with  that  total  impulse  would 
have  about  250  pounds  of  propellant  and  have  thrusts  ranging  from  6,000  to 
36,000  pounds.  The  highest  thrust  of  these  rockets  is  nearly  two  orders  of 
magnitude  too  low  for  the  PEA  application.  One  of  the  few  service  rockets 
which  have  a  thrust  near  the  1.44  million  pounds  required  is  the  TITAN  ill  C, 
which  is  ten  feet  in  diameter  and  85  feet  long--$omewhat  larger  than  desired. 
To  bridge  this  gap,  the  propellant  must  be  divided  into  many  small  pieces  so  ■ 
that  there  will  be  sufficient  burning  surface.  Some  very  strong  means  must  be 
devised  to  held  the  propellant  in  the  combustion  chamber  while  three  tons  per 
second  of  dense  gas  are  flowing  from  the  combustion  chamber  to  and  through  the 
nozzle.  The  drag  on  the  propellant  grains  from  this  enormous  flow  would  rip 
the  grains  loose  and  out  the  nozzle  before  they  were  burned.  The  highly 
divided  grains  and  the  distributed  How  area  required  would  make  such  a  rocket 
motor  much  larger  than  other  57,600  pound-secona  rocket  motors.  The  develop¬ 
ment  and  production  costs  for  such  rocket  motors  would  be  too  high  for  this 
program. 

The  problems  associated  with  the  two  high  thrust  rocket  motors  derive 
from  the  combination  of .  many  small  grains  of  propellant,  such  as  gun  propel¬ 
lant,  with  a  high  flow  through  the  bed  of  propellant  grains  and  out  a  nozzle. 
The  water  jet  PEA  alio  requires  highly  divided,  gun  like  propellant,  but  the 
combustion  chamber  is  not  vented.  The  burning  propellant  i$  always  enclosed 
in  a  rapidly  growing  volume  defined  by  the  anchor'  piston,  the  water  and 
launcher  tubes,,  and  the  gas-watet  interface.  Therefore,  the  propellant  grains 
do  not  have  to  be  restrained.  AH  that  is  needed  is  a  frangible  waterproof' 
container  that  can  keep  the  propellant  dry  during  the  deployment.  The  large 
scale  swaging  equipment  needed  to  make  the  water  tube  and  launch  tube  are 
available,  and  production  costs  are  moderate. 

■  The  conclusion  of  this  study  is  that  the  water  jet  PEA  is  clearly  the 
best  choice.  Experimental  development  of  this  system  should  be  started  as 
soon  as  possible. 


N»  S  t=  DISTANCE  BURNED 

'***.-—— >*"'  ODI=ODo-2* 

ID  =  ID  +  2  z 
H/-.H/-2. 

SURFACE  PER  GRAIN  «  it !OD„  +  ID  )<00_  -10  /2*-H  -4*1 

WHILE  *<(00  -ID  1/4 
o  o 

FIGURE  5.  THE  SURFACE  OF  A  MONOPERFORATED  GRAIN  AS  IT  BURNSAWAV 


1400 
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FIGURE  10.  WATER  VELOCITY  AT  EXIT  AND  AT  GAS  WATER  INTERFACE  AS  A  FUNCTION  OF  TIME. 
(PRESSURE  IS  ASSUMED  TO  BE  CONSTANT  AT  20,000  PSI) 


SY  (KPSI) 


FIGURE  14.  PLOT  OF  WATER  JET  TOTAL  SYSTEM  WEIGHT  VS  PRESSURE  AND  YIELD  STRESS,  SY. 

.  THIS  FAMILY  OF  CURVES  IS  BASED  ON  A  WATER  TUBE  INSIDE  DIAME  I  fcR  OF  1.25  FEET 


4.000 


FIGURE  15.  PLOT  OF  WATER  JET  SYSTEM  WEIGHT  VS  PRESSURE  AND  YIELD  STRESS,  SY.  THIS  FAMILY 
OF  CURVES  IS  BASED  ON  A  WATER  TUBE  INSIDE  DIAMETER  OF  1.5  FEET 


SY  <Kf>SI) 


(SCJNnOdl  XH9I3MW3XSAS 
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FIGURE  16.  PLOT  OF  WATER  JET  SYSTEM  WEIGHT  VS  PRESSURE  AND  YIELD  STRESS,  SY.  THIS  FAMILY 
OF  CURVES  IS  BASED  ON  A  WATER  TUBE  INSIDE  DIAMETER  OF  1.75  FEET 


(SCJNnOd)  1H0I9M  W3XSAS 


SY  (KPSO 


FIGURE  1».  PLOT  OF  WATER  JET  SYSTEM  WEIGHT  VS  PRESSURE  AND  YIELD  STRESS.  SY.  THIS  FAMILY 
OF  CURVES  IS  BASED  ON  A  WATER  TUBE  INSIDE  DIAMETER  OF  2.25  FEET 


ID-25 


FIGURE  19  PLOT  OF  WATER  JET  SYSTEM  WEIGHT  VS  PRESSURE  AND  YIELD  STRESS,  SY.  THIS  FAMILY 
OF  CURVES  IS  BASED  ON  A  WATER  TUtfE  INSIDE  DIAMETER  OF  2.5  FEET 


ID  «  2.75 


FIGURE  20.  PLOT  OF  WATER  JET  SYSTEM  WEIGHT  VS  PRESSURE  AND  YIELD  STRESS,  SY.  THIS  FAMILY 
OF  CURVES  IS  BASED  ON  A  WATER  TUBE  INSIDE  DIAMETER  OF  2.75  FEET 


.  7.000 


FIGURE  22.  PLOT  .OF  TOTAL  WATER  JET  PEA  WEIGHT  VS  WATER  TUBE  INSIDE  DIAMETER  AND  YIELD 
'•TRESS.  THIS  FAMILY  OF  CURVES  IS  BASED-ON  AN  OPERATING  PRESSURE  OF  25,000  PSI 


20,000  PSI 


47 


000  PSl 


5,000  PSI 
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APPENDIX  A 

RECOILLESS  COMPUTER  PROGRAM, 

INCLUDING  ITS  DESCRIPTION  AND  SAMPLE  CALCULATIONS 

The  computer  program  included  in  this  appendix  simulates  the  performance 
of  a  recoil  less  propellant  emolaced  anchor.  The  program  numerically  inte¬ 
grates  Equations  8  - 1 1 -5  through  14.  Equation  A-I-5  expresses  the  termination 
conditions,  arid  Equation  B-II-3  is  used  to  calculate  the  propellant  burning 
surface.  The  ideal  gas  equation  was  used  to  calculate  the  pressure  instead  of 
the  modifed  van  der  Waals  equation  (B-II-1),  because  the  value  for  the  covol¬ 
ume,  b,  was  not  known. 

The  main  program,  named  TETHER,  sets  values  for  the  constants,  and  reads 
in  a  file  (TAPE8)  of  operating  parameters.  Any  of  these  parameters  can  be 
changed  at  running  time  by  using  NAMELIST  STUFF.  After  the  input  parameters 
have  been  established,  the  main  program  calculates  auxiliary  parameters  and 
boundary  conditions.  All  of  the  input  and  most  of  the  calculated  parameters 
are  printed  so  that  the  user  will  have  the  the  input  parameters  and  the 
results  in  one  output.  The  main  program  then  prints  the  column  headings  and 
transfers  control  to  the  integration  subroutine.  With  the  aid  of  several 
other  subroutines,  it  integrates  the  equations  from  the  starting  conditions  to 
the  termination  condition  and  prints  the  results.  Control  is  then  returned  to 
the  main  program,  and  the  operator  is  queried  to  continue  or  stop.  If  the 
operator  elects  to  continue,  then'  he  is  requested  to  input  any  changed  data  ' 
via  the  NAMELIST  STUFF  statement.  The  entire  process  is  then  repeated  with 
the  new  set  of  input  data. 

•  The  integration  subroutine,  named  ROlNT,  can  integrate  any  number  of  1 
simultaneous , first  order,  di fferential  equations.  It  applies  the  Runge-Kutta 
method  for  numerical  integration.  From  the  main  program,  it  gets  sevtral 
parameters  such  as  the  integration  step  value  for  the  independent  varable, 
the  print  frequency,  the  initial . values,  and  the  termination  conditions.  For 
each  integration  step,  this  subroutine  must  transfer  to  the  OERIV  subroutine 
to  get  current  values  for  the  derivatives.  When  the  new  values  have  been 
calculated,  the  time  is  incremented  by  the  time  step,  and  the  new  variable 
values  are  stored  as  old  values.  This  process  Is  continued  unti 1  termination 
condition  is  met.  Not  every  step  is  printed.  The  first  and  last  ste>s  are 
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printed,  and  every  NPRINT  steps,  a  step  is  printed.  NPRINT  was  given  a  value, 
of  ten  in  this  program.  The  individual-  integration  steps  are  printed  by  a 
short  subroutine  named  PRINT. 

It  is  the  subroutine,  DERIV,  which  evaluates  the  equations  listed  in  the 
first  pararaph  of  this  appendix.  These  are  the  derivatives  which  are  inte¬ 
grated  in  RUINT.  The  pressure  is  also  evaluated  in  DERIV  by  use  of  the  equa¬ 
tion  of  state.  The  propellant  burning  equations  require  instantaneous  values 
of  the  total  area  of  the  burning  surface  of  all  the  propellant  grains.  This 
calculation  is  done  by  the  subroutine,  SURF.  The  version  of  SURF  listed  in 
this  appendix  calculates  the  burning  surface  of  a  number,  N,  of  seven  perf 
propellant  grains.  The  number  and  the  dimensions  of  the  grains  are  trans¬ 
ferred  to  the  function  subprogram  through  COMMON  statements.  The  instant¬ 
aneous  values  of  the  burned  distances  are  transferred  via  the  argument  list  to 
SURF.  The  derivation  of  the  equations  is  given  partly  in  Equations  B-II-3. 
This  gives  the  surface  of  the  grains  until  the  first  burnout,  when  the  seven 
perfs  have  grown  enough  to  become  tangent.  Subprogram  SURF  has  the  equations 
required  to  calculate  the  grains  to  complete  burnout.  While  these  equations 
aren't  very  complicated,  their  derivations  are  complicated.  Therefore,  the 
derivations  have  not  been  included. 

A  sample  data  set  and  a  computer  run  are  included  along  with  the  listing 
of  the  program.  These  show  the  way  the  data  must  be  arranged  and  the  way  that 
the  results  are  organized. 
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PROGRAM  TETHER ( INPUT, OUTPUT .TAPES ) 

REAL  MU.N.OD.ID.L.HA.MAV.ML.MLV.  JAY, STROKE 
LOGICAL  BEDONE 

DIMENSION  Y(30),NTERM(5),TERM(5) 

COMMON  /STUFF/RHOP.BEE.EN.TF.GAM.MU.N.OD.ID.L.PA.MA.FRIC.MAV, 
♦ML. MLV. ACS. APIS. CDD. AT, EPS, PRN.R. CD, CFVAC. STROKE 
EXTERNAL  DERI V 

DATA  GEE, JAY, R/32. 174, 773., 4970?./ 

NAMELIST  /STUFF/N.OD. ID.L,PINIT.V0,RH0P,BEE,EN,TF,GAf1,MU,D,MA, 
+FRIC, MAV, ML. MLV, ACS. APIS. CDD.AT.EPS. STROKE 
RE AD (8, 1 00 )N,OD. ID.L .PINIT , V0.RH0P.EEE.EN,TF,GAf1,hU.D,MA,FRIC, 
♦MAV, ML. MLV, ACS. APIS. CDD, AT, EPS, STROKE 
100  FORMAT (5F1 6. 8) 

15.0  READ  STUFF 

PA=D*63. 93+21. 16.  S 
G1 =(GAN-1 . I/GAM 

GFUNC=GA«/2.».3  ■  . 

GFUNC=GFUNC*‘(  -GFUNC/ < GAM- 1 . ) ) 

PRN= .  2 

200  PRNCLD-FRN 

PF  =  SQRT  (2.7(  GAM-1 .  )*U.-PRN**G1 )') 

PRN=(GFUNC/EPS/PF)**GAH 

IF ( AES ( ( FRN-PRNOLD )/PRN).GT.1.E-4)  GO  TO  200 
CFVA,C=EFS*  <2./G1*PRN**(1. /GAM ) -PRNMGAM+1 .)/ (GAM-1.)) 
AC=S3RT(GAM»R*TF/MU) 

CD=GAH*GFUNC/AC 

PRINT  300. TF ,GAH.MU,PRN,CFVAC,AC ,CB*GEE,RH0P*.01862.EN. 
♦B£E*12*144.**EN.N.0D*12.,ID*12.,L*12.,D.MA,M'AV,ACS. 


♦STROKE, VO. AT. APIS. CDD. ML, MLV, EPS, FRIC 


PRINT  350 

300  FORMAT (//"PROPELLANT  THERMODYNAMIC  DATA '/5X 'FLAME  TEMP  *' 

+  F3.2 "  DEG  R  V5X  'SPECIFIC  HEAT  RATIO  =  F6.4/5X  'MOLECULAR  UEIGHT  =  ' 

+F6.2/5X'EXIT  TO  CHAMBER  PRESSURE  RATIO  ='F6.5/5X 

♦''VACUUM  THRUST  COEFFICIENT  ='F6.4/5X'S0UND  SPEED  ='F10.2'  FT/S' 

+  /5X 'DISCHARGE  COEFFICIENT  ='F8.6'  LBM/S-LBF' 


♦/5X' PROFELLANT  DENSITY  ='F8.5 

♦'  LBH./CU  IN'/5X  BURNING  RATE  EXPONENT  -'F6.5/5X 

♦  'BURNING  RATE  COEFFICIENT  = ' F 1 0 . 8 '  IN/S-PSIA'//'GEOMETRIC  ' 

♦  PROPELLANT  DATA '/5X  'NUMBER  OF  GRAINS  ='F8.0/5X'GD  s'F/i.4'  IN' 

+  /5X  '  ID  *'F6.4'  IN'/jX' LENGTH  ='F8.5'  IN'// 

♦  ANCHOR  AND  LAUNCHER  D  AT  A ' /5X '  DEPTH  «'F3.0'  FTV5X 

♦  'ANCHOR  MASS  =  T4.2'  SLUGS VSX'ANCHOR  VIRTUAL  MASS  ~F6.2'  SLUGS' 

+  /5X 'ANCHOR  CROSS  SECTIONAL  AREA  ='F8 . 4'  SO  FT' 

+/5X'LAUNCHER  STROKE  =^8.4'  FT'/5X' INITIAL  FREE  VOLUME  *'F4.4'  FT' 
+  /5X 'ROCKET  THROAT  AREA  ='F8.5'  SO  FT'/5X'DRIVE  PISTON  AREA  *' 


♦  F8.r  SG  FT  '/5X  'ANCHOR  DRAG  COEFFICIENT  ='F8.A/5X 
♦'LAUNCHER  MASS  ='F6.2'  SLUGS'/SX 'LAUNCHER  VIRTUAL  MASS' =' 
♦F 6.2/  SLUGS'/SX'NOZZLE  AREA  RATIO  ='F8.4/5X 

♦  'FRICTION  FORCE  s'F10.2'  LBF './// ) 


NSWCTR  84-214 


350  FORMAT (30X  '< . VELOCITY - ; 

*  '  '  —  DISPLACEMENT  —  >  ' /oX  'TIME  UEB  BURN  GAS  MASS 

♦'ANCHOR  LAUNCHER  ANCHOR  LAUNCHER  GAS  VOL  GAS  TEMP  PRES' 


♦  •SURE  THRUST  '/5X'  (SEC)  '6X'  < IN) 

+  6X  (FT) ' 6X  (FT)  <CU  FT)  ( BEG  R: 
BO  400  1=1.10 
400  Y(  I )  =0. 

Y  C  3) =PINIT^VO ♦ HU/R/TF 


(SLUG) 
i  F  3 1 A  * 


<  r-t  r\  <  «  .  •  f  • 

Irvij  !•  l‘  o 


Y (3) =V0 
Y  ( ? )  =  T  " 


Y ( 1 3 ) =PINIT 

Y (161  =P’INI7  »AT*CF VAC  EPS  <- A T * P A 
NTERM (1 )  =  1 0 
TER.K  1 )=STROKE  • 

CALL  RUINT (1C,.  0001. 10..T..Y.1.N TERN. TERM. BERIV) 
REAB>, BEGONE  ■■  . 

IF ( .N.BEBONE)  GO  TO  150 

STOP 

END 


50 

100 

150 


200 

250 

300 


SUBROUTINE  RUINT(NU.'1.H.NPRINT.INIT,YHEU.NTC. NTERfl. TER.i.BERIV) 
dimension  yneu< 30 ). told ; :o,.b : so. 2< i :.g;;. so miterm; 30). 

♦  TERN (30) 

INTEGER  COUNT 
LOGICAL  IN  IT 

EQUIVALENCE! Y0LD( 1 ) .TOLD) 

DATA  B( 1 ) ,B( 2)/.5, .5/.D<3)  ,B( 4)/1 . .  J ./. COUNT. N/O.O/.B! 1  )/1 ./ 
Hri  =  -HH  , 


H  3  A  V  E  -  0 . . 

N  =  COUNT ♦ MOD (N-COUNT.NPR  INI) 

MTC  =  1 

IF!. NOT. INIT)  GO  TO  100 
CALL  PRINX(NUM.YNEU) 

COUNT =0 

N=NPRINT  , 

'hh=h 

DO  150  J=1,NUM 
YOLD(J)=YNEU( J). 

DO  250  J  =  1 ,4 

CALL  DERIV(NUM. YNEU,D,HH*B( J) ) 

DO  200  :=1 
G(  j.I)=DC)+HH 
YNEU 1 1 )  =  YULD ( I ) *0! J, 1 ) <B  ( J ) 

CONTINUE 
DO  3C0  1=1. HUM 

YNEU(  I  )=YCLD(  I  )KG(1, 1)10(4,1))  '6. ♦< 0 (2, 1  >  ♦(K 3 , 1 )  >/3. 
IF (MTC .GT.NTC )  GO  TO  380 
DO  325  J=MTC,NTC 
LTC=IA8S(NTERM( J) ) 

IF ( YNEU(LTC) *NTERM( J) .LT. 7ERH( J) ►NTERNT J) )  GO  TO  325 
MTC= JM  * 


4 
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HH=(TEKH(J)  -  YNE* ( L7C  5  5/0(110 
ri  5  A  V  E  -  H  3  A  U  E  *  H  H 
GO  T  G  103 
CONTINUE 
CONTINUE 

TF(HK.GT.O.)  GO  TO  400 
CAl.l  PRINT  (NUM. YNEUI 
:i!i=H3A'.'E 
TIT  URN 

"0" •• •  '•".•ii '  t  • 

o  L  u  i'  i  ‘  -  u  '1  I 

•  i  •  f'  »  * i .»•■•«*  t  n  cr,/> 

i  I  Li  J  J  •*  •  •  <(w  •  '  Ut'  *  ‘J  J  J 

"ALL  PRINT  i  N  U  fl .  Y  N  E  ii ) 

N=N*NPRINT 

GO  TO  50 


SUBROUTINE  DERI V t NUil . Y . B . DEL T ) 

REAL  il  U .  M  A ,  M  A  V ,  M  L .  M  L  V  ,J  A  V 
DIMENSION  Y(305.2( 30) 

COMMON  /STUFF •RHCr.SEE.EN.Tr .GAM.MU.N.CD.ID.L.PA.MA.rRIC.MAy. 
v'KL.  KLV.  ACS,  AF  IS.  CDD.AT.EPS.  PR’S.  R.CC.CrVAC.  STROKE 

jMI-'/Oi  '  ,  ,  , 

Y v  1 5 V=Y < 3 )  4R*Y(?)/mi/'Y(3) 

Y {  U )  =  Y  (1 5  HAT  *  { CFVAC  -1? 3  t-PA/Y  (13)) 

D( 2 )=EEE^Y( 1 5 ) ♦ »EN 
GP?-Rh'CPiSUEF(Y(2>>.4D(25 
IF  ■  UFR .EG. 0. )  D  (2 )  -  0. 

GEX  =  T(  1 5 ).»AT *L‘D'*SURT( TF/Yv? ) ) 

0*35 -G? R-GEX 


)  T  45  - ( APIS* i Y( ! 5  5  ~?A5  -  FRIC- .774  4  *AC3  *Y  ( 4)  f  » 2  *  CD  0  4  /  <  flA  frfl  AV ) 

)(3)-( APISH Yi  IS) -FA;  -FRIC-Y  ( 1 55  *AT*(CFVAC-Ei:'S *PA/Y< 15) ) ) / (ML *511.7  5 

!  '  •  i  -  V  1  '  - 


D  f  1  0 :  -  v  ( 4  i 

I  (3)=APIS»  (Y(4.5  i-Y(S5  )  ^PR/RHCP 

D ! ? )  5 ( Y ( 3 ) *Y ( 7 ! *G?R*CEL T *  TF ) / ( Y < 3 )  +GFR *  CELT  5  - (GAM -4 .  )+Y (?) / 
4 Y (354(3(3) -GPR/RHOP  5 


RETURN 

END 


SUBROUTINE  FRIN’-'N.Y: . 

CiMtNSiCN  i  ( -3 0  1 

PRINT  1 0 0 , f ( 1 ‘ . Y i 2 5 i 1 2 . . Y i 3 5 , Y ( 4 5 , Y ( 3 5 , Y ( o 5 , Y ( 2 ) . Y ( 3 5 , Y ( 7 ) , 
♦  Y ( 1 3 ) / T  4 4 . . Y ( 1  i 

rcR«At(3Fio.i.5Fi;.4.:ri'j.:.Fic'.(5) 
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FUNCTION  SURF (X 5 
KEAL  N.L, IB.GD.LO 

COMMON  /STUFF/RHGP.BEE.EN.TF.GAM.MU.N.OD. ID.LO.PA.MA.FRIC.MAV, 
+ML, MLV. ACS. APIS. CDD, AT, EPS, PRN.R, CD. CFUAC. STROKE 
P I  =  3  - 1 415926536 
RIO-IB/2. 

R00=0D/2. 

R I  =R 1 0  t-X 
RO-ROO-X' 

L=L0-,2.»X 
5= (RI+R01/2. 

SURF=0. 

IF (L.LE.O . .OR.RG/RI .LE.3./(5. -2.* SORT (3. ) ) '  RETURN 
SURF  =  2.*PI*(L*(R0  +  7.*RI)  +  R0**2-7.*RI*-*2)  *» 

IF CRI .LT.5/2. >  RETURN 
BETA=ASIN((5,*R0-3.*RI)/4./R0) 

GAMMA-BETA-PI/3. 

EPSILON=ACOS( CRO+RI )/4./RI ) 

IOTA=ASlN<RO/RI*COS (BETA) ) 

ALPHA=2.*PI/3. -EPSILON-IOTA 
SIGMA=PI/3.-2.*EPSILGN 

SURF =6. *(2*L*(RO*GAMMA+RI;*ALPHA)*RO*SIN( GAMMA >*SQRT ( (2.*RI*SIN( 
+  ALFHA/2. ) ) **2- ( RO*SIN( GAMMA) ) *«2)-RI**2*( ALPHA-SIN (ALPHA) )  ' 

♦♦.5*R0**2*<2.+GANNA-SIN<2.*GAMNA>>)*N 
IF(RI.GT.S/SQRT(3.))  RETURN 

SURF=SURF+<1.5*SQRT<3.>*S**2-?.*S*SaRT(RI**2-<S/2.>**2>  ' 
+-9.*RI**2*SIGMA+18.*RI*L*5IGHA)*N 
RETURN 
END 


TAPE3 


150.  ‘ 

,31666667 

.01666667 

.5 

864000 

1. 

3.23 

.0000044 

.8! 

5750 

1.22 

28.12 

300. 

200.  ■ 

1000 

50.. 

700. 

50. 

3. 

.1 

EOt  ENCOUNTERED. 

•  .5 

5.5' 

6. 

PROPELLANT  THERMODYNAMIC  DATA 
FLAME  TEMP  =  5750.00  DEG  R 
SPECIFIC  HEAT  RATIO  =1.2200 
MOLECULAR  WEIGHT  =  28.12 
EXIT  TO  CHAMBER  PRESSURE  RATIO  =.02678 
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APPENDIX  B 

DIRECT  ROCKET  COMPUTER  PROGRAM, 

INCLUDING  ITS  DESCRIPTION  AND  A  SAMPLE  CALCULATION 

While  the  mechanism  of  the  direct  rocket  PEA  is  significantly  different 
from  that  of  the  recoilless  PEA,  their  performances  dan  be  simulated  by  very 
similar  computer  programs.  The  recoil  less  PEA  has  two  parts  which  move 
differently,  although  interdependent^.  The  anchor  moves  forward  under  the 
force  of  the  gas  pressure  acting  on  a  piston,  less  the  internal  and  external 
frictions.  The  recoil  less  launcher  is  accelerated  backward  by  the  same  force, 
a  pressure  on  a  piston.  This  force  is  diminished  by  the  same  small  internal 
friction  that  affects  the  anchor  and  by  the  very  large  rocket  thrust.  The  net 
force  acting  on  the  launcher  should  be  vey  small  as  compared  to  the  force 
acting  on  the  anchor.  Since  the  net  force  acting  on  the  launcher  is  a  small 
difference  between  two  large  forces,  that  small  net  force  could  be  positive  or 
negative.  A  positive  force  would  cause  the  launcher  to  accelerate  to  the 
rear;  a  negative  force  would  accelerate  it  forward.  Regardless  of  the  sign  of 
the  force,  its  magnitude  will  be  reduced  by  the  external  friction.  The  effec¬ 
tive  masses  of  the  anchor  and  the  launcher  are  both  increased  by  their  virtual 
masses.  The  mass  of  the  launcher  is  steadily  decreasing  because  of  the  pro¬ 
pellant  gases  being  ejected.  The  temperature  of  the  gases  inside  the  combus¬ 
tion  chamber  and  piston  must  be  computed  continuously  because  the  gas  is  doing 
expansion  work  and  new  propellant  gas  is  continuously  being  mixed  with  the 
• Ider  gas. 

The  direct  rocket  moves  as  a, single  unit;,  therefore,  only  one  velocity 
and  displacement' must  be  integrated.  There  is  no  internal  drag.  The  force 
accelerating  the  whole  system  is  equal  to  the  rocket  thrust  less  the  external 
hydrodynamic  force.  The  net  mass  is  increased  by  the  virtual  mass  and 
decreased  by  mass  of  the  propellant  gas  exhausted.  The  propellant  temperature 
remains  constant  inside  the  combustion  chamber,  because  the  gas  does  no  expan¬ 
sion  work  until  it  reaches  the  nozzle. 

Because  the  details  of  the  calculations  are  so  similar  to  those  of  the 
recoilless  PEA,  it  was  most  efficient  to  modify  the  recoil  less  program  to 
handle  the  direct  rocket.  The  SURF  subprogram  was  changed  to  handle  single 
per?  propellant  grains  because  \he  direct  rocket  does  not  need  a  progressive 
burn. 


,?  * 
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PROGRAM  DIRSOC (INPUT, OUTPUT. TAFES) 

REAL  MU .  N .  OB .  I B .  L .  M  A .  MAV .  JAY  .STROKE 
LOGICAL  BEDGNE 

DIMENSION  Y»30).MTERH(5).TERH(5) 

COMMON  /STUFF /SHOP . BEE ,EN ,  TF  ,G  Ail  ,HU.  N ,  OD .  IB ,  L .  PA  ,KA, MAV, 

♦ACS  .ODD.  AT,  EPS,  PRii .  R .  CD ,  CFV  AC.  STROKE 
EXTERNAL  DERIV 

DATA  GEE. JAY. R/32. 174, 773. .4970?./ 

NAMELIST  /STUFF /N . OD . I D . L . PIN! T . VO .RHQP , DEE , EN . TF . GAM . HU . E , HA . 

♦ MAV. ACS. ODD. AT. EPS. STROKE 

READ(3,100)N.QB,IB.L,PINIT.VCsRH0?.BEE.EN.TF,3AM,MU.D,i1A. 

.  +MAV, ACS . CBB . AT . EPS. STROKE 
TOO  FORMAT ( 3F t 6 . 3 ) 

150  READ  STUFF 
PI=4.'*ATAN(1 . ) 

RA*B*43. 98+21 14.8 
G1 -(GAM-1 . )/GAM 
GFUNC-GAM/2.+.5 

GFUNC=GFUNC*'*(-GFUNC/(GAN-1.5) 

PRN=,2 

200  PRNQLD=PRN 

PF=SQRT(2. /(GAM-1. L*(1.-PRN**G1)> 

PRN=<GFUNC/EPS/PF>**GAM 

IF(A8S((PRN-PRNCLD)/PRN).GT.1.E-4)  GO  TO  200 
CFVAC=EPS*(2./G1*PRN**(1 . /GAM ) -PRN*( GAM+1 .)/ (GAM-1 .)) 
AC=S3RT(GAM*R*TF/HU) 

CD=GAH*GF'UNC/AC 

PROP.M=N*RHOP*P I  *L* ( CD**2--ID**2 )  /4 . 

PR  I  NT  300.TF.GAM.I1U.  PRN ,  CFVAC .  AC .  CD  i  GEE .  SHOP* .  0 1 342 .  EN . 

+  BEE+.1 2*1  44.**EN,N.QB*12.,ID*12..L*12.  .D.MA.MAV.PRQPM.ACS  . 
♦STROKE. 90, AT.CDB.EPS 
'  PRINT  350 

300  FORMAT  (//"PROPELLANT  THERMODYNAMIC  BATA '/5X 'FLAME  TEMP  -' 

+F8.2'  DEG  R'/3X' SPECIFIC  HEAT  RATIO  * 'f6 . 4 /3X ' MOLECULAR  HEIGHT  * 
♦F4.2/3X'EXIT  TO  CHAMBER  PRESSURE  RATIO  =  'F4.S/SX' 

♦  'VACUUM  THRUST  COEFFICIENT  *'F4.4/3X'S0U*B  SPEED  * 'FT 0.2'  FT/S' 

♦  /Sa'DIS.CHARGE  COEFFICIENT  *'F8.4'  LBM/S-LBF ' 

♦  /SX 'PROPELLANT.  DENSITY  *'FS. S 

+  '  LBM/CU  IN '/SX  'BURNING  RATE  EXPONENT  t'Fft’.S/SX 
♦'BURNING  RATE  COEFFICIENT  »'F10.8'  JN/$~PSIA'//'GE0MET8IC  ' 
♦'PROPELLANT  BATA '/5X 'NUMBER  OF  MONO? ERF  GRAINS  •'F8.0  • 

: ♦/SX'DB  •*>F4.4'  IH'/SX'IB  *'F6.4'  I N'/3X 'LENGTH  *'FS.5'‘  IN'// 

♦  'SYSTEM  BATA-/3X 'DEPTH  *'F5.0'  FT'/SX 

♦'SYSTEM  MASS  *'f4.2'  SLUGS '/3X' SYSTEM  VIRTUAL  MASS  ♦'F4.2'  SLUGS 
♦/SX'PROPELLANT  MASS  *'F30. 4'  SLUGS' 

♦/SX 'SYSTEM  CROSS  SECTIONAL  AREA  *'F3.4'  S3  FT' 

♦ /S X 'STROKE  *'FS.4"  FT'/SX 'INITIAL  FREE  VOLUME  *'F9.4'  FT' 

• ♦/SX'ROCKET  THROAT  AREA  «'F8,5'  $8  ft'‘  ' 

♦/SX'ANCKCS  DRAG  COEFFICIENT  *'F8.4 
♦/SX'MOEELS  AREA  RATI3. *'F3.4///> 
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35C  FORMAT  (  oX  'TIME  USB  BURN  GAS  HASS  FRCP  HASS' 

+  ■"  VELOCITY  DISTANCE  GAS  VOL  PRES,' 

* ' SURE  THRUST ' /5X'{SEC)'6X' ( IN)  (SLUG)  (SLUG) 
■  »oX'(FT)  (Cli  ri)  (PSIA)  {BOUNDS)'//) 


l  r i/ll 


*?  w  V  i  *  *  !  ~  ‘V  • 

Y  \  3 ) =  P I  »•*  I T  *  V  0  >  f*  m  /  X  /  T  F 


Y <  1  o) i;ti  i  -Aw  , 

ntzr.i:d=6 

TERHi I )=STROKE 

CnLii  Nu i N T  ( /  •  •  o 00 1  ■  iO«  •  t  • « i  . » *iMuAHf  i  Ei\H « i>Ei\  t  v ) 

PRINT  500 

500  FORMAT i // 'DO  YOU  UANT  TO  STOP  {.!.)' OR  CONTINUE  (.7.1  COMPUTING? ' ) 


READ  ’>■ .  EEOCNE 
Ir  ;.N.EE52NE)  30  TO  ISO 


SUBROUTINE  mUINT.'NUH.H.NPPINT.INIT.YNEU.NTC.NTERH.TERH.DERIU) 
0 1  .TENSION  YNt'«  ( 30  - .  YOLO  ( 30 ) .  D 30  5 . 8.(4 Q  ( 4 . 30  ) . NTERH  ( 30 ) . 

+  TERM (30) 

INTEGER  COUNT 
LOGICAL  IN  IT 

H  «  *1  *•  \  I  1  •  J  *.  1  '  ,  J  •  ,  Z  ■  •  6  \  3  •  •  3  l  ^  ’  *  •  •  *  •  /  ,  C0L’?«T  .  .v  /0  •  0  /  •  D  C  i  %i  /  1  »/ 
IiH-  -V.' 


*?  *  -  ■>  v  .*  ♦  *  v  i.*  (•«**■  v.  C  o  N  T  •  t*  r\\  I N  T  / 

r.:c^ 

!*:  r 


M’Vl  Jv  -V 

w  n  *  r  •*  1 N  *  *  !l  'J  ^ f  ^  t  *  . 


•  •  *  . 


•  **  v  «.  ^  t  *  a  « i  v  *  •  Y  »*  t  ■*  •  J  *  •  •  *  i ' 


.•  W  ..  y  . 


00  Y‘iC.:  i;  ..VOL o:  I  ■  0  ’ I  •  .0; 
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MTC-J4- ■ 

HH* i TZSJl ( J )  -  Y.1EU C LTC ) » / B i  LTC  > 


",  ;«  .  i.-|i  JM‘t,  '  f 

r  ft  »n  « 

O  U  I  v>  IVV 

1_-  u»'^  «  irtlJL 

J  f*  I  ill'wU 

*  f  »  «1  M  **■  ▼  !*■  \ 

x  r-  \  ■  O  i  •  V  . 

*  1  !  or 

DLL  I  I 

S  Z  T  U  S  * 


IF  ? C C li N I •  N E •  K  i  u ’j  •  -  ■*'' 
CALL  PRINT  (SUM.  THK) 
N-NH*  PRINT 


iu 


c  rt  if 


SUBROUTINE  DfRIVvNU.I.Y .D.DELT  > 

REAL  M«.,1A.“A'V .  JAY 
DIMENSION  Y(30>.D(?0> 

COMMON  /STUFF/RNCP.BEE .EN. TF. GAM, MU. N. 01 
f ACS. cur. AT. EPS. FRN.R. CD, CFVAC. STROKE 

i  »  v  -  V»9 
WM  t  */  r  w  • 

T { 135  -  Y<3)  TF 

Y ( 1 4 1  =  r  < 1 3 5 1  A T  M  CT VAC  EPS  1  PA/ Y  !15)) 
D(2)=BEE»  UT35  **EN 
GPRsRrOP*S'JRF (Y>‘2) )  0(0 
IF (3PK.Z3.3.  00--0. 

3EX  -  Y  >'  1')  ‘AT  iCD 
»(3)sGPR-5cX 


DUU-GtX 

drag=.??»uacsu;5>*  i2.cr:2 
D(35  =  (  U15)  *  AT  * :  CFVAC  •  ZTZ-  ’<  PA/Y  ( 
D(6)=Y<3) 

D(7)=G?R‘RHCP 

RETURN 

END 


SUBROUTINE  PRINT vN.Y) 

DIMENSION  YCO- 

PRINT  lOC.Yd  i.YCJ)  1 1L-.. US ).Y(  15. 

iyuo u  ;a.,y(:s) 

FORMAT ( 3F 1 0. i.  IF  1 0 . 4.F-1 0.2 ,F1 0 .0 ' 

RETURN 

END'  ' 


c r-  C  J  «'E  Cu  Cd  Cc:  in.  tc!  ...J  CO  *•«  CO  rc: 
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DO  YOU  WANT  TO  STOP  (.T.)  OR  CONTINUE  (.F.)  COMPUTING? 
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APPENDIX  C 

WATER  JET  COMPUTER  PROGRAM, 

INCLUDING  ITS  DESCRIPTION  AND  A  SAMPLE  CALCULATION 

The  flow  equations  and  their  derivatives  were  not  included  in  the  main 
text  because  their  derivation  is  too  lengthy  and  complex.  This  appendix  also 
includes  the  computer  program  which  evaluates  these  equations  and  gives  a 
detailed  view  of  the  events  that  happen  during  the  .040  seconds  of  a  launch. 
The  input  data  and  a  typical  run  are  included. 

Only  two  of  the  Navier  Stokes  Equations  are  needed  because  the  flow  is 
incompressible.  The  two  equations  needed  are  the  equations  of  conservation  of 
mass  and  momentum.  The  equation  for  the-  conservation  of  mass  is 

C-l  p  A(z)  U(z,t)  *  p  A(x)  U(x,tl  . 

Since  the  flow  is  incompressible,  the  densities  on  both  sides  of  Eqn.  C-l  are 
equal  and  can  be  canceled.  It  is  more  convenient  to  express  the  flow  areas, 
A(z)  and  A(x),  as  the  product  of  the  main  barrel  flow  area,  AQ,  and  a  function 
G(x)  or  G(z).  The  term,  G(z),  gives  the  flow  area  at  any  point, along  the 
length,  z,  as  a  fraction  of  AQ.  With  the  water  tube  shape  discussed  in  the 
text,  G ( z )  has  a  value  of  1  through  the  straight  portion  of  the  tube  and  less 
than  one  in  the  convergent  nozzle.  G(x)  represents  the  area  at  x,  the  posi¬ 
tion  of  th£  gas-water  interface.  Again,  both  sides  have  equal  terms,  AQ,  and 
they  can  be  canceled.  This  simplifies  Eqn.  C-l  to 

C-2  G{z)  U(z,t)  =  G(x)  U(x,t)  . 

The  one-dimensional  momentum  equation, is 

C-3  i  +  u(z,t)  ■  0  . 

Rearranging  Eqn.  C-2  gives 

C-4  U(z,t)  -  U(x,t)  . 
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Taking  a  partial  derivative  with  respect  to  z  gives 

C-b  ^  =  -  U(x,t)  G(x)  , 

3z  G  (z) 

where  G'(z)  =  d^z-^-  . 


Taking  a  partial  derivative  of  U(z)  with  respect  to  t  gives 


The  reader  should  be  warned  that  U(x,t)  is  not  standard  notation.  It  is  used 
to  signify  the  water  velocity  at  any  time  where  z=x.  That  explains  why  there 
is  no  derivative  with  repect  to  z  and  why  the  derivative  with  respect  to  time 
is  not  a  partial  derivative.  Substituting  Equations  C-4,  5,  and  6  into  Equa¬ 
tion  C-3  gives 


c_7  I  aP-Ckt!  +  G(xj  +  u2,x  G ‘  (x)_ 

C  '  p  3Z  dt  G(z )  U  'X*t>  GU) 

-  U2(x,t)  G2(x)  =  0  . 

Gd(z) 

Now  multiply  Eqn.  C-7  by  dz  and  integrate  over  z  from  x  to  L. 


♦  U2(x,t )  G2(x) 


*  *!»>  ^  JT  b& 

Lg2{L)  G2 ( x ) J 


+  u  (x,t )  G'(x) 


Now  solve  Eqn.  C-B  for  .  , 


'  dU(x,t)  . 
dt 


P(x,t)  -  P(l,t)  -  U2(x,t)  |l  -  V1- 


P  G{x) 


‘  dz 

5T?T 


Eqn.  C-9  can  be  integrated  with,  respect  to  time  to  give  the  velocity  of 
the  gas-water  interface  as  a  function  of  time.  This  velocity,  in  turn,  can  be 
integrated  with  respect  to  time  to  give  the  position  of  the  gas -water  inter¬ 
face  as  a  function  of  time.  In  order  to  integrate  Eqn.  C-9,  one  needs  to  be 
able  to  calculate  the  pressure  P  at  x  and  l  for  all  times  in  the  domain  of  the 
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calculation.  Similarly,  one  needs  a  function  for  G(z)  so  that  G(x),  G(L), 
G'(x),  and  the  integral  of  the  reciprocal  of  G(z)  from  x  to  L.  G(z)  has  a 
constant  value  of  one  (1)  for  values  of  z  from  zero  to  L$;  therefore,  G'(z) 
equals  zero  over  that  range  of  z.  The  value  of  G(z)  decreases  monotonically 
from  one  to  the  minimum  value  as  z  goes  from  L$  to  L.  Consequently,  the  value 
of  G‘(z)  is  negative  in  this  range.  Because  G(z)  is  determined  by  two  differ¬ 
ent  functions  in  the  two  regions,  it  is  necessary  to  evaluate  the  integral  of 
the  reciprocal  of  G(z)  in  two  steps.  (This  breakup  should  only  be  used  when  x 
is  less  than  L$.) 

c'10  JxgT5T  '  £S W)  '  "here  x<Ls- 

G(z)  has  a  constant  value  of  one  (1)  while  z'  is  in  the  interval  from  zero  to 
L$.  Therefore,  the  first  integral  on  the  right-hand  side  of  Eqn.  C-10 
simplifies  to  the  integral  of  dz  over  the  same  range.  That  integral  reduces 

where  x  v  l.  . 

When  x  is  greater  than  L$,  the  integral  is  not  broken  up.  It  is  evaluated  by 
integrating  directly  from  x  to  L. 

To  calculate  x,  the  distance  traveled  by  the  gas-water  interface,  it  is 
necessary  to  integrate  the  following  equations. 

C - 12  .  ^  =  U(x,t) 

The  boundary  conditions  for  the  velocity  and  the  distance  are 

U(x,t)  *  0  when  t  *  U 

x  *  0  when' t  *  0 

As  with  the  other  two  PEA  options,  it  is  r.ecessary  to  be  able  to  calcu¬ 
late  the  pressures,  P(x,t)  and  P(L,t).  In  an  all-up  analysis,  the  pressures 
would  be  calculated  from  first  principles  as  was  done  for  the  recoilless  and 
direct  rocket  PEAs.  This  was  not  done  for  the  water  jet  PEA  computer  pro¬ 
gram.  The  pressures  were  set  to  constant  values.  Even  though  the  actual 


C- 13 
C-14 
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pressure  would  not  be  precisely  constant,  it  would  not  be  very  far  off.  This 
simplification  saved  time  and  permitted  many  runs  to  be  compared  without  the 
secondary  effects  of  pressure  variation  blurring  the  relationships. 

There  is  another  shortcoming  in  the  analysis  as  it  is  presented  here.  No 
account  is  taken  of  the  mass  of  the  propellant  gas  in  the  dynamic  equations. 
Some  of  the  energy  from  the  expanding  gas  must  be  used  to  accelerate  the  y'-s 
as  well  as  the  water.  The  density  of  the  gas  is  about  one  tenth  that  o',  water 
when  the  pressure  is  10,000  psi  and  two  tenths  when  the  pressure  is  i:0,000 
psi.  The  added  mass  of  the  gas  is  insignificant  at  the  start  of  the  cycle 
when  most  of  the  tube  is  filled  with  water.  Although  the  quantity  of  gas  is 
greater  during  the  main  mid-portion  of  the  run,  it  is  also  insignificant  there 
because  the  acceleration  is  zero  then.  The  only  time  that  it  has  a  signifi¬ 
cant  effect  is  very  near  the  end  of  the  cycle,  when  there  is  very  little  water 
and  the  mass  of  the  gas  exceeds  that  of  the  water.  Only  the  last  part  of  the 
cycle  is  significantly  affected  by  the  mass  of  the  gas.  The  acceleration  of 
the  interface  would  be  very  high  and  approach  infinity  as  the  denominator  of 
the  main  term  in  Equation  C-9  goes  to  zero;  In  reality,  the  mass  of  the  gas 
would  prevent  this  from  happening,  but  that  mass  does  not  show  in  Equation  C- 
9.  The  computer  program  in  this  appendix  handles  this  in  a  nonrigorous,  but 
effecti  ve,  way. 

Equation  C-9  gives  the  acceleration  of  the  gas  water  interface.  If  one 
multiplies  the  numerator  and  denominator  of  the  main  term  by  the  term  Ag,  then 
one  can  see  that  that  first  term  is  a  force  divided  by  a  mass  (an  accelera¬ 
tion).  In  the  enclosed  program,  the  mass  of  the  gas  is  given  as  the  product 
of  the  mass  density  of  the  water  and  the  effective  length  of  the  gas  column 
(Again  AQ  has  been'  factored  out.)  In  the  program,  this  mass  of  gas  pen  unit 
area  is  added  to  the  mass  of  water  per  unit  area  in  the  denominator  of  the 
main  term. 
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PROGRAM  PEA (INPUT, OUTPUT) 

COMMON  L.LN.GL, RHC .  PX . PL , AO 
DIMENSION  Y(30)  .NTERrli  5)  .TERM (5) 

REAL  I.L.LS.LN.NINTS 
LOGICAL  REDONE 
EXTERNAL  DERIV 

NAMELIST  /AG A  IN/PX , DEPTH, L.LN .GL . AO .TSTEP.NPRINT 

RHO= 1 . 7545 

BEOQNEs.F. 

PRINT  50 

50  FORMAT ( '1  TYPE  IN  VALUES  FOR  PX. DEPTH, L.LN.GL .AO, TSTEP.NPRINT') 
READ  t.PX. DEPTH, L.LN.GL, AO, TSTEP.NPRINT 
ICO  LS=L-LM 

PL=14.7*DEPTH*.4347 
DO  200  J-I.30 
200  Y ( J )  - 0 . 

PR IN :  300 .PX. PL . DEPTH , L , LN . AO ,GL*AO . TSTEP 
■  300  FORMAT (  '  ‘  INPUT  DATA '  /SX  'DRIVE  PRESSURE  ='F10.2'  P3IA  V5X 
A '  AMBIENT  PRESSURE  *  'F10.2"  PSIAV5X  DEPTH  =  'F6.0'  FEET  V5X 
B  TUBE  LENGTH  ='F3.4'  FEET  '/5X' NOZZLE  LENGTH'FS. 4'  FEET  V3X 
C  'CHANNEL  CfiCSS  SECTION  *'F8.4'  SQ  FEET  '/5X  'NOZZLE  AREA  «'F8.4 
D'  S3  FEET " / 5 X  '  INTEGRATION  STEP  =  'El  0.3'  SEC 
E///6X  TIME  DISTANCE  '6X'V(X) '6X'V(Y )  THRUST  IMPULSE'/ 

F5X  C3ZC)  (FEET)  (FT/S)  (FT/S)  (POUND)  (LB-S) '//> 
NTERM(n=3 

TERNM  < -.'>?? H  ,  ' 

CALL  RUINT(4, TSTEP.NPRINT ,. T .. Y . I .NTERM, TERM. DERIV) 

PRINT  400 

400  FORMAT (.'////)  '  .  , 

REAP  AGAIN 

IF ( . NOT . BEDOi i)  GO  TO  100 
END 


SUBROUTINE  RU I N  T ( NUM , H . NPR I NT , IN I T , YNE  U , N  T  C , NTERM , TERM .DERIV) 
DIMENSION  YNEU ( 30 ) . TOLD( 30) , D( 30) ,B( 4 ) ,Q( 4 , 30 ) . NTERM (30) , 

♦  TERM( 30 ) 

INTEGER  COUNT  ~  '  '  ■ 

LOGICAL  I MI T 

EQUIVALENCE ( YOlDX 1 ) , TOLD) 

DATA  B(U,B(2)/,5*. 5/,B(3),B(4)/1.,1 ./, COUNT, N/0,0/,D< 1 )/1 ,/ 
HSAVE=0. 

N  =  COUHT  +  MOD(N- COUNT.  NPRINT) 

M  T  C  =  1 

IF ( .NOT. INIT )  CO  TO  100 
CALL  PRINT (MUM, YNEU) 

COUNTS 
M= NF  R I NT 
MO  HH*»! 

100  DO' ISO  J* 1 .MUM 


"  C-5 
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150  Y0LD(J)=YNEU(J> 

DO  250  J  =  1 ,4 

CALL  DFRIV(NUM, YNEU.D) 

DO  200  1=1, NUM 
Q( J, I ) =D( I ) +HH 

200  YNEU ( I) = YOLD ( I ) +0 ( J . I ) *B ( J ) 

250  CONTINUE 

DO ' 300  1=1, NUM 

300  YNEU(  I  )=Y0LD(  I )  +  (Q( 1 , I  )+Q(4, 1 ) )/6.  +  (0(2. 1 ) *G ( 3, 1 <  5/3. 
IF(MTC.GT.NTC)  GO  TO  380 
DO  375  J=MTC,NTC 
LTC=IABS(NTERM( J) ) 

IF ( YNEU(LTC)*NTERM( J) .LT.TERM( J)*NTERM( J) )  GO  TO  375 
HTC=J+1 

HH=(TERH( J)-YNEU(LTC) >/B(LTC) 

HSAVE=HSAVE+HH 
GO  TO  '00 
375  CONTINUE 

380  CONTINUE 

,IF(HH.GT.O.)  GO  TO  400 
CALL  PRINT  (NUM. YNEU 5 
HH=HSAVE 
RETURN 

•100  COUNT  sCOUNT  +  1 

IF  (COUNT. NE.NJGO  TO  50 
CALL  PRINT  (NUM.  YNEU) 

N=N+NPRINT 
GO  TO  50 
END 


SUBROUTINE  KRIWINUH. Y.D) 

COMMON  L,LN, GL.RHO.PX.F'L. AO 
REAL  L.LN.LS 
DIMENSION  Y (30) .0(30) 

D  ( 1 )  =  1 . 

CALI  N02 ( Y< 3 )  ,GPR IME, GINT  .RAT ) 

D(2).=<(PX-PL>*144./RH0+.5*Y(2)**2*<1  ♦  *($X/6t)**2-R'Af*U  .*SX**2) )) 
A/GX/GINT-Y(2)**2*GPRIME/GX 
D(3)=Y(2) 

D( 4)=RH0*A0*GX*< (L*Y( 3) )*D(2)*Y(2)+*2* (GX/Gl  -  I . > ) 

Y(20)=D<4) 

Y(21)=!Y(2)*GX/5L 

RETURN 

END 


K>  * 


NSWC  TR  84-214 


SUBROUTINE  NQKX.GX  ,GPRIME,GINT,RAT) 

CGMMON  L,LN.GL.RHO,PXtFL,AO 

REAL  L.LN.LS 

LS=L-LN 

IF (X.GT.L*. 999999999)  X=L* .999999999 
PI=3. 1 4159 
RAT=PX*6.24E~6 
GPRIilE  =0 .  • 

GX  =  1 . 

•  G I N  T  =  R  A  T *X+LS -XH.N/SQRT ( GL ) 

IF(X.LE.LS)  RETURN 

GX= <  1 .  f GL ) /2 . ♦  ( 1 . -GL ) n . *COS ( P I* ( X-I.S )  /LN ) 

GPR I H£= ( GL - 1 . )*PI/2./LN*SIN(PI*(X-LS)/LN) 

TERN=2.*LN/PI/SQRT(GL)*ATAN(SGRT(6L>*TAN(PI*(X-LS)/2./LN) 

GINT -RAT  KLS+TERH)  +  LN/ SORT (GL ) -TERM 

RETURN 

END 

SUBROUTINE  PRINT (NUM.Y) 

DIMENSION  Y(30) 

PRINT  '100, Yd)  ,Y(3),Y(2),r (21). Y(20),T(«) 

100  F0RHAT(F1 0.6.F1 0. 4.2F1 0.3.2E10 .3) 

RETURN 

END 


GO 

TiPE  IN  VALUES  FOR  PX, DEPTH, L.LN.GL. AO, TSTEP.NPRINT 
0000 . , 300 . . 7 . , 1 . , - 1 . 2 . 65 . . 0000 1 , 1 00 
INPUT  DATA 

DRIVE  PRESSURE  *.  20000,00  PSIA 
AMBIENT  PRESSURE  =  145.71  PSIA 

DEPTH  =  300.  FEET 

TUBE  LENGTH  7.0000  FEET 
NOZZLE  LENGTH  1.0000  FEET 
CHANNEL  CROSS  SECTION 
■  NOZZLE  AREA  *  ,2650  SO  FEET 

INTEGRATION  STEP  ■*  .100E-04  SEC 


2.6500  SO  FEET 
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TIME 

DISTANCE 

V  (X ) 

(SEC) 

(FEET) 

(ET/S) 

0 .  CP  ■'(.'MO 

0.0000 

0.000 

.001000 

.0705 

125.667 

.002000 

.2203 

163.966 

.003000 

.3886 

170.709 

.004000 

.5600 

171.733 

..005000 

.7318 

171.881 

.006000 

.9037 

171.902 

.007000 

1  .0756 

171.904 

.008000 

1.2475 

171 .905 

.009000 

1.4194 

171.905 

.010000 

1.5913 

171 .905 

.011000 

1.7632 

171.905 

.012000  ' 

1.9351 

171.905 

.013000 

2.1070 

171.905 

.014000 

2.2739 

171.905 

.015000 

2.4508 

171.905 

.016000 

2.6227 

171.905 

.017000 

'2.7946 

171.905 

.018000 

2.9665 

171.905 

.019000 

3.1384 

171.905 

.020000 

3.3104 

171.905 

.021000 

3.4823 

171.905 

.022000 

.3.6542  ’ 

171.905 

.023000 

3.3261 

171.905 

.024000 

3.9980 

171 .905 

.025000 

4.1699 

171.905 

,026000 

4.3418 

171.905 

.027000 

4.5137 

171.905 

.028000 

4.6856 

1 71 .905 

.029000 

4.8575 

171.905 

.030000 

5.0294 

17,1 .905 

.031000  1 

5.2013 

171.905 

.032003 

5.3732  ' 

171,905 

.033000 

5.5451 

171 .905 

.034000 

5,7170 

•171.903 

.035000 

5.8389 

171.905 

.036000 

6.0610 

173.317 

.037000 

6,2431 

,196.648 

.033000 

6. 4769 

296.886 

.038303 

6.9930 

2235.200 

vm 

THRUST 

IMPULSE 

(FT/3) 

(FOUND) 

(LB-S) 

0.000 

0. 

0.  ' 

1256.668 

.342E+07 

•434E+04 

1639.658 

. 1 77E+07 

.725E+04 

1707.092 

. 1 44E+07 

.881E+04 

1717.334 

. 1 39E+07 

. 1 02E+05 

1718.811 

. 1 38E+07 

. 1 1 6E+05 

1719.016 

.13BE+07 

. 1 30E+05 

1719.044 

. 1 38E+07 

. 1 43E+05 

1719.047 

.138E+0? 

.  157E+05 

1719.048 

.  1 38E+07 

. 1 71 E+OS 

1719.048 

. 1 33E+07 

.1 85E+05 

1719.0*8 

.138E+G7 

.  1 99E  *05 

1719.048 

.  1 38E+07 

.21 2E+05 

1719.048 

. 138E+07 

.226E+05 

1719.048 

.  1 38E+07 

.240E+05 

1719.043' 

.  1 33E+07 

.  254IE  4-05 

1719.048 

.  1 33E+07 

.247E+05 

1719.048 

.  1 38E+07 

.281E405 

1719.048 

.  1 38E+07 

.295E+05 

1719.048 

.  133E+07 

. 309E+05 

1719.048 

.  1 38E+07 

.323E+05 

1719.048 

.  1 38E+07 

.336E+05 

1719.048 

.  1 33E+07 

.350E+05 

1719.048 

. 1 33E+07 

•344E+05 

1719.048 

. 1 38E+07 

.378E+05 

1719.046 

. 1 38E+07 

,3?1E*05 

1719.048 

, 1 38E+07 

.405E+05 

171.9.048 

.138EI07 

.417005 

1719.043 

.  1 38E+07 

.433E+05 

1719.048 

.13SE+07 

.447E+05 

1719.048 

,138E<-07 

.  460E+05 

1719.048 

.138E+0? 

.474005 

1719.048 

.13SE107 

,483E*05 

1719.048 

.133E+07 

.302EM55 

1719.048 

.138E+07 

•515E+05 

1719.048 

.133E+07 

.529E+05 

1719.094 

.142E+0? 

.S43E»05 

1720.661 

,  131E+07 

.558E+05 

1729.629 

.163E+07 

•373E+05 

2257.625 

.428E+03 

.536E+05 

NSWC  TR  84-214 


APPENDIX  D 

SYSTEM  SIZE  AND  WEIGHT  COMPUTER  PROGRAM, 

INCLUDING  ITS  DESCRIPTION  AND  A  SAMPLE  CALCULATION 

The  BASIC  computer  program  included  in  this  appendix  calculates  a  matrix 
of  weights  and  dimensions  of  water  jet  PEA  launchers.  The  program  was  coded 
to  assume  that  all  of  the  water  jet  launchers  had  to  deliver  5. 76  x  ID4  foot¬ 
pounds  of  impulse  at  a  thrufc  level  of  1.44  x  10^  pounds.  It  also  assumes 
that  all  of  the  launchers  are  made  of  steel  with  a  weight  density  of  492 
pounds  per  cubic  foot.  Within  these  constraints,  the  weights  and  dimensions 
of  the  water  jet  PEA  launcher  were  calculated  for  a  matrix  of  values  of  the 
inside  diameter  of  the  water  tube,  the  driving  pressure,  and  the  yield  stress 
of  the  structural  material. 

Four  equations  from  the  main  text  are  needed  to  calculate  the  dependent 
variables.  Equation  E-I-l  is  needed  to  calculate  the  inside  diameter  of  the 
launch  tube;  Equation  E-I-b  to  calculate  the  weight  of  the  launch  tube;  Equa» 
tion  E -I -10  to  calculate  the  length  of  the  water  tube;  and  Equation  E-I-13  to 
calculate  the  weight  of  the  water  tube.  The  total  system  weight  calculated  in 
this  program  does  not  include  the  weights  of  the  anchor,  piston,  propellant, 
or  water.  It  is  the  sum  of  the  weights  of  the  launch  tube  and  the  water  tube. 
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00100  MARGIN  132 
0 0 1 10  R1  =  t.9837 
001  20  READ  13, T.  RO 
00130  READ  PI ,P2,P3 
00140  READ  D1.D2.53 
00150  READ  SI  ,S2,';3 


00160  P9=4  »A  TN  n  > 


0017,0  FOR  D-D1  TO  D2  STEP  53 
00180  A=P9*D*2/4 

00190  PRINT  “INSIDE  DIAMETER  3F  UATCR  TUBE  ="P"  FEET. 

00200  PRINT  "WATER  TUBE  CROSS  SECTIONAL  AREA  SQUARE  FEET* 
00210  PRINT  "TOTAL  IMPULSE  =“I3"  POUND  SECONDS." 


00230  PRINT 

00240  PRINT  TAB (21 ) ;"UATER  TUBE" 

00250  PRINT  "  P-f'A  DARREL  ID  IEN3TK 

00260  PRINT  "■>- -r--- . - . . TOTAL  S'?3TEi  UEI"; 

00270  PRINT  "rHT .  " 

00280  PRINT  “  TPS I )  (FEET  •  <PEET'“; 

00290  PRINT  *  .  -rPOUNOS";  . 

00300  PRINT  -) . •• 

00310  PRINT 

00320  PRINT  -HELP  STRESSES  'PSD . 

00330  FOR  S*SI  TO  S2  STEP-  S3 
00340  PRINT  USING  ,  S/1  4  »; 

00350  NEXT  S 
00360  PRINT 
00370  PRINT 

00380  FOR  P-P1  TO  P2  STEP  P3 


00390  V=I3/S0R'2»R1  MP-T/AM 


00400  L*V/A*D/3 

0(410  PRINT  USING  00420.  P/1  *  IJSPR 'T 1 1.T/P° » ;l ; 

00420  :  lllllltll 

00430  FOR  S-SI  TO  S2  STEP  S3 


00440  U0*R0»12»I/(3-2*P) 

00450  U=R0»’P‘DM2»'l  •  D  )/'*'$ -2  * f- >*  0/ ( S -1 .5*P>>  »U0 

00460  PRINT  USING  00470, UJ 

00470  :  mull 

00480  NEXT  5 

00490  PRINT 
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00500  NEXT  P 
00510  PRINT 

00520  PRINT  T AB ( 21 ) ;"UATER  TUBE" 


00530  PRINT  H  P-PA  BARREL  IB  LENGTH  M; 

00540  PRINT  H<— . — - - - - BARREL  UE!" 

00550  PRINT  "GHT . . . . >" 

00560  PRINT  "  (PSI)  (FEET)  (FEET)”: 

00570  PRINT  M  '  < . - - - ' - (POUNDS 

00580  PRINT  M> . . . . - . . >" 

00590  PRINT 


00600  PRINT  "TIELD  STRESSES  (PSI) - >"J 

00610  FOR  S=S1  TO  S2  STEP  S3 
00620  PRINT  USING  ”##333###" ,S/1 44; 

00630  NEXT  S 
00640  PRINT 
00650  PRINT 

00660  FOR  P -PI  TO  P2  STEP  P3 
00670  V=I3/SQR(2*R1KP-T/A)) 

00680  L=V/A+D/3 

00690  PRINT  USING  00700,  P/144;S0R(T*4/P/P?) JLJ 
00700  :#########  #«###.###  ######.## 

00710  FOR  S=S1  TO  S2  STEP  S3 
00720  UO-RO* 1 2 * T  '(S~2*P) 

00730  PRINT  USING  00740, UO; 

00740  :  ####### 

00750  NEXT  S 
00760  PRINT 

00770  NEXT  P  ,  ■ 

00780  PRINT 
00790  PRINT 
00800  NEXT.  B 

00810  BATA  57600,1440000,492 
00820  BATA  1 440000 ,3600000 , 1 44000 
00830  BATA  1.25, 3, .25 
00840  BATA  1.44  E7,2.88  E7.1.44  E6 
00850  END 


D-3 


INSICE  OIANETER  OF  U^TEK  USE  ‘  2.25  FEET.  WATER  TUBE  CROSS  SECTIONAL  AREA  =  3.97603  SQUARE  FEET 
TOTAL  IMPULSE  =  S^iOO  POUND  SECONDS. 

IHRUST  =  1440000  FOUNDS. 


4183  3604  3166  2323 


HELD  STRESSES  <* SI' . :--1  100000  HOOOO  120000  1  30000  1*0000  150000  160000  170000  180000  190000  20*000 
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